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1  Abstract 

The  thorough  understanding  of  the  formation  and  the  relaxation  of  the  plasma  produced 
in  the  shock  layer  developed  during  the  re-entry  of  a  spacecraft  in  the  upper  layers  of  the 
earth’s  atmosphere  is  crucial  in  order  to  prevent  damaging  of  its  outer  surface.  Among 
the  different  points  to  be  studied,  the  chemical  aspects  are  particularly  important:  the 
mechanical  characteristic  time  scale  of  the  flow  being  short,  the  flow  is  indeed  in  chemical 
nonequilibrium.  In  addition,  the  inner  storage  of  energy  of  the  different  species  of  the 
flow  leads  to  other  mechanical  behaviors  than  those  observed  in  classical  low  temperature 
flows.  These  nonequilibrium  effects  have  to  be  taken  into  account  accurately  in  codes 
devoted  to  the  study  of  the  interaction  between  the  plasma  and  the  surface. 

In  this  paper,  we  propose  to  focus  our  attention  on  the  detailed  kinetics  of  an  air  plasma 
under  nonequilibrium  conditions.  First,  we  present  a  time-dependent  collisional-radiative 
(CR)  model  for  atomic  oxygen  to  study  in  detail  ionization  and  three-body  recombination 
rates  of  oxygen  in  high  temperature  air  plasma  flows.  Second  we  present  a  time-dependent 
CR  model  for  air  taking  into  account  13  species  and  numerous  excited  states  and  working 
over  a  wide  range  of  pressure  and  temperature.  For  typical  conditions  encountered  in  re¬ 
entry  flows,  we  compare  this  detailed  CR  scheme  for  air  with  simplified  schemes  (proposed 
by  Park,  Dunn  &  Kang  and  Gupta  et  al.)  usually  implemented  in  high  enthalpy  air  flow 
codes.  Finally,  a  ID  code  simulating  the  stagnation  point  boundary  layer  near  the  wall  of 
the  spacecraft  is  presented  and  the  influence  of  the  gas  chemistry  on  the  wall  catalycity 
is  discussed. 


2  Introduction 

During  the  re-entry  of  a  spacecraft  into  the  upper  layers  of  the  atmosphere  of  a  planet,  the 
flow  is  hypersonic  and  a  shock  layer  is  formed  in  front  of  the  vehicle  (Anderson  (1989)). 
In  crossing  the  shock  layer,  the  flow  is  hardly  slowed  down  and  the  gas  is  strongly  heated 
(up  to  15000  A'  (Carlson  and  Hassan  (1992)))  and  therefore  becomes  a  plasma.  Near  the 
surface,  a  boundary  layer  is  formed  due  to  the  difference  between  the  temperature  inside 
the  shock  layer  and  the  wall  temperature  and  this  induces  high  energy  flux  densities  to  the 
vehicle  surface.  These  fluxes  have  to  be  reduced  to  prevent  damaging  of  the  spacecraft. 
To  choose  correctly  the  material  of  the  thermal  shield  it  is  necessary  to  study  the  thermal 
interaction  between  the  gas  medium  and  the  vehicle  (Hankey  (1994))  and  therefore  to 
study  thoroughly  the  plasma  formed  in  the  shock  layer  (i.e.  different  species  concentra¬ 
tions  close  to  the  wall,  their  internal  degree  of  excitation,  etc).  In  the  present  context, 
the  flow  is  far  from  equilibrium;  after  an  ionization  phase,  the  plasma  is  recombining  as  it 
approaches  the  thermal  shield  and  the  velocity  is  high.  As  a  result,  the  chemistry  is  time- 
dependent  in  the  frame  of  the  moving  plasma.  Its  calculation  needs  an  accurate  kinetic 
scheme.  Moreover,  as  the  pressure  varies  over  a  wide  range  during  the  re-entry  phase, 
the  kinetic  scheme  has  to  take  into  account  a  sufficient  number  of  binary  and  ternary 
collisions.  Currently,  despite  the  large  number  of  experimental  and  theoretical  studies  on 
high  temperature  air  chemistry,  uncertainties  remain  on  the  accuracy  of  numerous  reac¬ 
tion  rates. 

In  particular,  the  effective  ionization  and  recombination  rates  of  atoms  are  difficult  to 
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determine  experimentally  and,  then  most  available  results  have  been  derived  from  theo¬ 
retical  studies.  However,  a  reaction  such  as  0+  +  2e_  <-»■  O  +  e~  results  from  numerous 
different  elementary  processes  between  the  atomic  levels  of  oxygen.  Therefore,  the  ac¬ 
curacy  of  theoretical  effective  ionization  and  recombination  rate  coefficients  depends  on 
the  choice  of  the  atomic  model  and  on  the  choice  of  the  rate  coefficients  for  the  elemen¬ 
tary  processes.  In  section  3,  we  present  a  time-dependent  collisional  radiative  model  for 
atomic  oxygen  in  order  to  determine  effective  ionization  and  recombination  rate  coeffi¬ 
cients  for  1018m~3  <  ne  <  1021m"3  and  8000/1  <  Te  <  20  000/1 .  We  will  show  that  a 
time-dependent  approach  allows  to  follow  the  relaxation  of  atomic  level  populations  and 
to  study  the  validity  of  the  quasi-steady-state  hypothesis. 

As  far  as  we  know,  no  time-dependent  collisional-radiative  model  dedicated  to  re-entry 
applications  exists  for  air,  taking  into  account  N2,  02,  NO,  N,  O,  N%,  ,  NO+ ,  N+ , 

0+,  02-,  O-  on  their  ground  state  and  on  numerous  electronic  excited  states.  Although 
some  collisional-radiative  models  have  been  earlier  elaborated  by  Teulet  et  ah  (2001)  and 
Sarrette  et  al.  (1995),  no  model  works  between  1  kPa  and  the  atmospheric  pressure.  In 
addition,  considerable  improvements  have  been  made  about  data  in  term  of  excitation  by 
electron  impact,  dissociative  recombination  of  NO+,  02  and  and  their  branching  frac¬ 
tions  as  well  as  vibrational  processes.  In  section  4,  we  propose  to  focus  on  two  important 
points  in  the  case  of  the  plasma  formed  during  the  re-entry  into  the  earth’s  atmosphere: 
(1)  the  time-dependent  chemistry  of  an  ionizing  or  recombining  air  plasma  and  (2)  the 
comparison  with  existing  simplified  kinetic  schemes  used  up  to  now  in  Navier-Stokes  codes. 
This  study  is  based  on  the  elaboration  of  a  time-dependent  collisional-radiative  model  for 
heavy  particles  temperature  T 4  and  electron  temperature  Te  less  than  10000  K.  The  re¬ 
sults  of  the  elaborated  CR  model  are  finally  discussed  for  different  re-entry  conditions. 
Finally,  in  section  5,  we  present  the  influence  of  the  choice  of  the  kinetic  scheme  for  air  on 
the  wall  heat  flux  on  a  TPS  material  in  the  stagnation  point  boundary  layer  for  different 
conditions. 


3  Collisional-radiative  model  for  O  atoms 

In  this  section,  we  propose  to  set  up  a  collisional  radiative  model  for  atomic  oxygen  in 
order  to  determine  effective  ionization  and  recombination  rate  coefficients  for  1018m-3  < 
ne  <  1021m~3  and  8000/1  <Te  <  20  000/1.  In  an  earlier  study,  Taylor  and  Ali  (1986)  de¬ 
termined  these  rate  coefficients  for  slightly  different  conditions:  1022m~3  <  ne  <  1025m~3 
and  leV  <  Te  <  3eV.  Other  previous  collisional  radiative  models  for  oxygen  have  been 
devoted  to  the  study  of  the  population  densities  at  the  stationary  state  (Gomes  et  al. 
(1990);  Soon  and  Kune  (1990))  and  relaxation  times  (Cacciatore  and  Capitelli  (1976)). 
Since  the  study  of  Taylor  and  Ali,  numerous  experiments  and  computations  have  been 
carried  out  to  improve  the  accuracy  of  atomic  cross  sections  and  Einstein  coefficients  in 
oxygen.  In  this  work,  we  have  tried  to  determine  the  best  possible  coefficients  available  in 
the  literature.  The  data  critically  compiled  by  Soon  and  Kune  (1990)  have  been  used  with 
a  few  updates  on  e— O  inelastic  collisions  (Doering  (1992))  and  on  Einstein  coefficients 
(Biemont  and  Zeippen  (1992);  Bhatia  and  Kastner  (1995)).  In  this  study,  we  consider 
a  spatially  uniform  and  electrically  neutral  plasma.  We  assume  that  the  distributions 
of  energies  of  particles  are  Maxwellian  and  that  the  atom-atom  and  ion-atom  inelastic 
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collisions  are  negligible.  The  treatment  of  radiation  is  limited  to  either  optically  thin  or 
optically  thick  cases. 

Unfortunately,  no  experimental  data  are  available  on  the  effective  ionization  and  recombi¬ 
nation  rates  of  atomic  oxygen  to  compare  with  our  results.  In  high  temperature  air  kinetic 
schemes  used  in  plasma  flow  codes,  different  analytical  expressions,  which  accuracy  is  un¬ 
known,  are  proposed  for  these  rate  coefficients  (Park  (1990);  Gupta  et  al.  (1990);  Losev 
et  al.  (1994)).  Therefore,  this  work  is  the  opportunity  to  study  the  validity  of  these  simple 
analytical  expressions  for  the  effective  ionization  and  recombination  rates  of  atomic  oxy¬ 
gen.  Furthermore,  the  time-dependent  approach  used  in  this  paper  allows  us  to  follow  the 
relaxation  of  atomic  level  populations  and  to  study  the  validity  of  the  quasi-steady-state 
hypothesis. 

3.1  Atomic  model 

In  this  work,  forty  coalesced  energy  levels  of  the  neutral  oxygen  atom  have  been  considered 
(Moore  (1993)).  Their  characteristics  (for  each  level  index:  level  configuration,  energy 
and  statistical  weight)  are  given  in  Table  1.  As  in  the  nitrogen  atom,  there  is  a  large 
energy  gap  in  the  oxygen  atom  between  the  two  first  metastable  levels  (2s22pA  4D  and 
2s22p4  4S)  and  the  other  excited  levels  which  lie  in  a  narrow  region  of  energy.  However, 
it  is  interesting  to  note  that  in  the  oxygen  atom,  the  resonant  level  is  only  the  fifth 
level,  located  at  9.52  eV  from  the  ground  state.  In  their  atomic  model  for  oxygen,  Soon 
and  Kune  (1990)  considered  only  the  nine  first  levels  of  our  model.  In  this  case,  the 
last  level  lies  1.69  eV  under  the  ionization  limit.  With  our  model,  this  energy  gap  is 
only  of  0.124  eV.  The  influence  of  the  atomic  model  on  the  results  will  be  discussed  in 
the  following  sections.  In  the  electron  temperature  range  of  our  study,  we  consider  only 
the  ground  ionic  term  of  0+,  2s22p3  4S  (13.62  eV)  which  is  taken  as  the  ’’core”  for  all 
atoms.  Processes  involving  0~  ions  could  a  priori  also  influence  the  determination  of  the 
ionization  and  recombination  rate  coefficients.  However,  Soon  and  Kune  (1991)  showed 
that  for  Te  >  8000/i ,  the  contribution  of  0~  ions  to  production  of  electrons,  positive  ions 
and  excited  atoms  is  negligible.  Therefore,  in  this  study,  the  influence  of  0~  ions  has 
been  neglected. 

3.2  Elementary  rate  coefficients 

In  a  uniform  plasma,  the  evolution  of  the  population  density  of  an  excited  atomic  level  i 
is  due  to  different  elementary  processes 


Oi  +  e  ^  Oj  +  e  for  i  <  j  [Cl3,  Cji] 


(1) 


where  CtJ  (m3  s~4)  and  Cji  (m3  s_1)  are,  respectively,  the  electron-impact  excitation  rate 
coefficient  for  the  transition  from  the  level  i  to  the  level  j  and  its  inverse  deexcitation  rate 
coefficient 


Oi  +  e  0+  +  2e 


(2) 


where  (3iC  (m3  s_1)  and  /3Ci  (m6  s”1)  are,  respectively,  the  electron  impact  ionization  rate 
coefficient  from  the  level  i  and  the  three-body  recombination  rate  coefficient  on  the  same 
level 


Oi  +  hv  for  i  <  j  [A,,;] 


(3) 
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i 

Is2  + 

Ei(e  V) 

9i 

1 

2s22p4(3P ) 

0.000 

9 

2 

2s22pA(lD) 

1.970 

5 

3 

2s22pA(1S) 

4.190 

1 

4 

2s22p33s(5S°) 

9.150 

5 

5 

2s22p33s(3S°) 

9.520 

3 

6 

2s22p33p(5P) 

10.740 

15 

7 

2s22p33p(3P) 

10.990 

9 

8 

2s22p34s(5S°) 

11.840 

5 

9 

2s22p34s{3S°) 

11.930 

3 

10 

2s22p33d(5D°) 

12.090 

25 

11 

2s22p33d(3D°) 

12.100 

15 

12 

2s22p34p(5P) 

12.300 

15 

13 

2s22p34p(3P) 

12.370 

9 

14 

2s22p33s'(3D°) 

12.550 

15 

15 

2s22p35s(5S°) 

12.670 

5 

16 

2s22p35s(3S°) 

12.710 

3 

17 

2s22p33s'(1D°) 

12.740 

5 

18 

2s22p34d(5D°) 

12.760 

25 

19 

2s22p34d(3D°) 

12.770 

15 

20 

2s22p34f(5F,  3F) 

12.780 

56 

21 

2s22p35p(5P) 

12.860 

15 

22 

2s22p3{3P) 

12.890 

9 

23 

2s22p36s(5S°) 

13.030 

5 

24 

2s22p36s{3  S°) 

13.050 

3 

25 

2s22p35d(5D°,  3D°) 

13.080 

40 

26 

2s22p35f(5F,  3F) 

13.087 

56 

27 

2s22p36p(5P) 

13.130 

15 

28 

2s22p36p(3P) 

13.140 

9 

29 

2s22p37s(5S°) 

13.220 

5 

30 

2s22p37  s(3  S°) 

13.230 

3 

31 

2s22p3[6d(5D°,  3D°) 

+6/(5F,  3F) 

+6 g{5G°,  3G°)] 

13.250 

168 

32 

2s22p38s(5S'°) 

13.330 

5 

33 

2s22p38s(3S'°) 

13.340 

3 

34 

2s22p3[7d(3D°,  5D°) 

+7d(5F,  3F)} 

13.353 

96 

35 

2s22p39s(3S°) 

13.412 

3 

36 

2s22p38d(5D°,  3D°) 

13.418 

40 

37 

2s22p310s(3S°) 

13.459 

3 

38 

2s22p39d(3D°) 

13.464 

15 

39 

2s22p3lls(3S°) 

13.493 

3 

40 

2s22p319d{3  D°) 

13.496 

15 

41  (0+  ion) 

2s22p3(AS°) 

13.620 

4 

Table  1:  Atomic  level  model. 


9-6 


RTO-EN-AVT-142 


Detailed  and  Simplified  Kinetic  Schemes  for  High 
Enthalpy  Air  Flows  and  Their  Influence  on  Catalycity  Studies 


where  (s  3)  is  the  transition  probability  (Einstein  coefficient)  from  the  level  j  to  the 
level  and 


where  afR  (m3  s^1)  and  afD  (m3  s^1)  are,  respectively,  the  radiative  and  dielectronic 
recombination  rate  coefficients  on  the  level  i.  Dielectronic  recombination  is  the  result  of 
a  two-stage  process 


(5) 


Oi  +  hu 


where  O*  represents  an  autoionizing  state  of  the  O  atom. 


3.2.1  Electron-impact  excitation  and  deexcitation  rate  coefficients 

Currently,  no  complete  set  of  precalculated  cross  sections  exists  for  oxygen.  Therefore, 
for  low  atomic  levels,  measurements  of  cross  sections  have  been  used,  when  available. 
For  the  1  — ►  j  <  7  transitions,  the  cross  sections  measured  by  Doering  et  al.  (  Doering 
and  Gulcicek  (1989b),  Doering  and  Gulcicek  (1989a),  Gulcicek  et  al.  (1988)  and  Gulcicek 
and  Doering  (1988))  have  been  used  in  agreement  with  the  discussion  of  Soon  and  Kune 
(1990).  For  the  1— >-2  transition,  it  is  interesting  to  note  that  the  shape  and  magnitude  of 
the  cross  section  have  been  confirmed  recently  by  Doering  (1992).  For  the  1— >  9<  j  < 
19  transitions,  the  theoretical  and  experimental  cross  sections  reviewed  by  Laher  and 
Gilmore  (1990)  have  been  used.  For  all  the  other  transitions,  excitation  cross  sections 
have  to  be  determined  theoretically.  In  this  work,  we  have  used  for  low  atomic  levels  the 
Vainshtein  formalism  which  is  certainly  more  adapted  than  hydrogenic  models  retained  in 
earlier  studies.  It  would  be  interesting  to  use  cross  sections  obtained  with  more  accurate 
quantum  mechanical  methods  such  as  close-coupling  and  /?- matrix  methods.  But  these 
methods  are  computer  intensive,  and  have  been  completed,  up  to  now,  only  for  a  few 
systems  (Pradhan  (1988)).  For  ’’difficult”  transitions  (i.e. ,  involving  two  atomic  electrons), 
a  binary-encounter  approximation  is  used.  All  details  are  given  in  the  paper  of  Soon  and 
Kune  (1990).  For  atomic  levels  close  to  the  continuum,  all  atoms  are  almost  hydrogenic. 
Therefore,  the  formalism  initially  developed  for  hydrogen  by  Gryzinski  (1965)  is  used  for 
these  levels,  and  is  no  doubt  more  adapted  than  the  Vainshtein  formalism. 

The  use  of  different  theoretical  models  may  result  in  abrupt  changes  in  the  elementary 
rate  coefficient  values.  In  order  to  eliminate  this  anomaly,  already  mentioned  in  other 
plasma  studies  (Park  (1971))  we  introduce,  in  the  same  way,  a  correction  factor  such 
that  elementary  rate  coefficient  values  vary  smoothly  from  one  model  to  the  other.  As 
an  example  Fig.  1  shows  the  rate  coefficient  for  the  1  —>  j  transition.  Reliable  rate 
coefficients  based  on  the  review  of  Laher  and  Gilmore  (1990)  are  available  up  to  j  =  19. 
Gryzinski  rate  coefficients  up  to  j  =  19  have  higher  values  but  we  note  that  in  both  cases, 
the  rate  coefficient  decreases  rapidly  up  to  j  =  10  and  decreases  more  slowly  for  higher 
values.  Therefore,  in  our  model,  rate  coefficients  for  1  — >  j  >  19  are  derived  using  the 
Gryzinski  rate  coefficients  scaled  using  the  rate  coefficient  of  the  1  — >  19  transition.  We 
have  checked  that  discrepancies  on  the  effective  ionization  and  recombination  rate  values 
due  to  this  correction  remain  within  a  factor  of  2. 
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Figure  1:  Rate  coefficients  of  the  1  — >  j  transitions.  Open  circles: 
theoretical  and  experimental  rate  coefficients  reviewed  by  Laher  and 
Gilmore  (1990).  Dashed  line:  Gryzinski  (1965)  rate  coefficients. 
Solid  line:  our  model. 


3.2.2  Electron-impact  ionization  and  three-body  recombination  rate  coeffi¬ 
cients 

For  the  ground  state  and  the  two  lowest  excited  levels,  the  models  retained  by  Soon  and 
Kune  (1990)  are  used.  For  the  other  atomic  levels,  we  have  used  the  classical  model 
approximation  derived  by  Gryzinski  and  Kune  (1986). 

3.2.3  Spontaneous  emission 

For  excited  levels  the  energy  of  which  is  less  than  or  equal  to  11.93  eV,  the  transitions  and 
the  Einstein  coefficients  reviewed  by  Soon  and  Kune  are  used.  However,  for  the  5— >1  and 
9— >1  transitions,  the  Einstein  coefficients  determined  recently  by,  respectively,  Bicmont 
and  Zeippen  (1992),  and  Bhatia  and  Kastner  (1995)  have  been  used.  For  transitions 
issuing  from  levels  lying  higher  than  11.93  eV,  the  best  possible  Einstein  coefficients 
available  in  the  literature  have  been  used  (Bhatia  and  Kastner  (1995);  Fuhr  and  Wiese 
(1990);  Tayal  and  Henry  (1989);  Doering  et  al.  (1985);  Verner  et  al.  (1996)). 

3.2.4  Radiative  and  dielectronic  recombinations 

The  analytical  expressions  proposed  by  Nussbaumer  and  Storey  (1983)  are  used  to  cal¬ 
culate  the  dielectronic  recombination  coefficients  for  the  effective  (direct  plus  cascade) 
transitions  to  the  terminating  levels  i  =  1  and  4.  The  total  dielectronic  recombination 
rate  is 

aRD  =  £  aRD  (6) 

*=1,4 

As  in  nitrogen,  radiative  recombination  coefficients  are  significant  only  for  the  effective 
transitions  to  the  three  first  atomic  levels.  Using  the  same  approach  as  Soon  and  Kune 
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(1990),  we  have  calculated  and  fitted  the  associated  total  radiative  recombination  rate, 
and  obtained 

3 


aRR  =  =  2.12  x  io-18T-°-29  m3s-1 

2—1 


(7) 


3.3  Master  equations 

Taking  into  account  the  different  elementary  processes  mentioned  above,  the  rate  equa¬ 
tions  for  the  production  of  an  excited  atom  on  the  level  i  can  be  written  as 


dNi 

dt 


=  Ni  =  V'  rieNjCji  +  T  NjAjiKji  +  neN+  [aRRKRR  +  aRD krd  +  nef3ci\ 


j>i 


-N, 


^  ^  neCij  T  ^  ^  +  ne/3j, 


k<i 


(8) 


for  1  <  i  <  np,  where  np  is  the  total  number  of  atomic  levels,  ne  (m-3)  the  electron 
number  density,  Nt  (m-3)  the  population  density  of  the  atomic  level  i,  and  N+  (nr-3)  the 
number  density  of  0+ .  For  the  escape  factors  denoted  Kij,  krr  and  krd  ,  respectively, 
for  bound-bound,  free-bound,  and  dielectronic  radiation,  we  consider  only  the  limit  cases 
where  they  are  equal  either  to  0  (optically  thick  plasma)  or  to  1  (optically  thin  plasma). 
The  rate  equation  for  the  production  of  electrons  is 

=  rie  =  Yl  nlN+Pci  -  neN+  [a^R^R  +  afD^RD]  (9) 

i  i  i 

Finally,  the  electrical  neutrality  of  the  plasma  imposes  ne  =  N+ .  We  denote  Na  —  ^  Ni, 
the  total  atomic  density  and  Nt  —  Na  +  N+,  the  total  particle  density.  The  system  of 
differential  equations  [Eqs.  (8)  and  (9)]  has  been  integrated  numerically  in  using  LSODE 
Hindmarsh  (1980).  In  each  set  of  calculations,  Te  and  Nt  are  constant. 


3.4  Results 

3.4.1  General  behaviour  of  time-dependent  population  densities 

Figure  2  shows  the  temporal  evolution  of  the  level  populations  for  a  typical  recombining 
situation  at  t  =  0  s:  ne  =  1022m~3,  Te  =  12  000A',  and  bi  =  Ni/N*  =  10-4,  where  N* 
is  the  Saha  density  of  level  i.  The  plasma  is  considered  to  be  optically  thin  and  all  the 
elementary  processes  described  above  are  considered.  This  example  clearly  puts  forward 
that  different  time  scales  characterize  the  global  relaxation  process.  For  t  <  10“ 10  s,  in 
this  example,  Ni,  N2,  IV3  and  ne  vary  only  slightly  in  comparison  to  other  excited  state 
population  densities.  For  5  x  10“9s  <  t  <  10~6s,  the  latter  reach  a  plateau  whereas  Ni, 
N2,  N3  and  ne  vary  significantly.  Finally,  all  the  system  tends  toward  a  stationary  state 
(reached  at  f  ~  10“3s  in  this  example).  It  is  interesting  to  note  that  the  populations  of 
excited  states  at  the  plateau  are  very  different  from  those  at  the  stationary  state.  It  is 
important  to  mention  that  the  general  characteristics  of  the  results  are  independent  of 
the  initial  ionizing  or  recombining  condition. 
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Figure  2:  Temporal  evolution  of  densities  for  an  initial  recombining  situation:  at  t  =  0  s, 
Te  =  12  000  K,  ne  =  1020nr~3  ,  and  6*  =  Ni/N*  =  10-4,  where  N*  is  the  Saha  density  of  level 
i.  All  elementary  processes  are  taken  into  account,  and  the  plasma  is  considered  to  be  optically 
thin.  Numbers  denote  atomic  level  numbers.  Atomic  level  population  densities  are  divided  by 
their  statistical  weights. 


Figure  3:  For  the  same  initial  recombining 
condition  as  in  Figure  2.  Solid  line:  \rie\. 
Solid  line  with  open  circles:  N±.  Dashed 
line:  N\  +  N2.  Dot-dashed  line:  N\  +  N2  + 
N3. 


Figure  4:  For  the  initial  recombining  con¬ 
dition  of  Fig.  2.  Solid  line:  |ne|/n3.  Solid 
line  with  open  circles:  IVi/n3.  Dashed  line: 
(Ni  +  A^/n3.  Dot-dashed  line:  (IV1  +  IV2  + 
lY3)/n3. 
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Figure  3  shows  the  temporal  evolutions  of  |rie|  and  Nt  for  1<  i  <3.  For  t  >  rmjn  (rm;n  ~ 
5x10”9s  in  this  example),  we  have 


Ni  ~  0  for  i  >  3  (10) 

and  then  a  quasi-steady-state  condition  of  the  form 

3 

\rie\  =  (11) 

i= 1 

can  be  written.  This  condition  has  been  already  established  by  Cacciatore  and  Capitclli 
(1976)  in  oxygen.  For  hydrogen  as  shown  by  Bates  et  al.  (1962),  the  corresponding  quasi- 
steady-state  condition  takes  into  account  only  the  ground  state  of  the  atom.  In  oxygen 
and  also  in  nitrogen  as  shown  by  Bourdon  and  Vervisch  (1996),  it  is  necessary  to  take  into 
account  the  contribution  of  the  two  first  low-lying  levels.  In  oxygen,  the  fourth  level  is 
metastable  too  and  then  one  could  have  expected  to  have  for  oxygen  a  quasi-steady-state 
(QSS)  condition  of  the  form 

4 

\rie  |  =  y (12) 

1=1 

In  the  electron  temperature  and  density  range  of  our  study,  we  have  checked  that  Eq. 
(11)  is  the  quasi-steady-state  condition  for  oxygen. 


3.4.2  General  method  to  determine  effective  ionization  and  recombination 
rate  coefficients 

For  the  global  reaction  0+  +  2e~  O  +  e~ ,  we  define  Ay,  the  effective  three-body 
recombination  rate  coefficient,  and  kt,  the  effective  ionization  rate  coefficient,  such  that 


dne 

dt 


=  nP 


killeNa  -  krnl 


(13) 


Both  rate  coefficients  have  to  be  time  independent,  but  may  depend  on  electron  temper¬ 
ature  and  densities.  To  determine  kr,  the  method  is  to  consider  at  t  =0  s  a  recombining 
situation.  In  this  case,  in  early  times  the  first  term  on  the  right  hand  side  of  Eq.  (13) 
is  far  much  smaller  than  the  second  one.  Then,  if  kr  exists,  there  is  a  plateau  in  the 
temporal  evolution  of  the  ratio  rie  \ /nj  and  kr  is  the  value  of  this  ratio  at  the  plateau. 
Conversely,  to  determine  A;*,  one  has  to  consider  at  t  =0  s  an  ionizing  condition  and  to 
study  the  temporal  evolution  of  the  ratio  rie/{neNa).  ki  is  the  value  of  this  ratio  at  the 
plateau.  For  the  initial  recombining  situation  of  the  last  section,  Fig.  4  shows  that  a  value 
of  Ay  can  be  determined  for  Tmin  <  t  <  Tmax ,  with  Tmax  ~  10~5s.  The  minimal  time  rmm 
corresponds  to  the  necessary  time  to  reach  a  quasi-steady-state  (Figs.  2  and  3)  and  is 
denoted  Tqss  in  the  following.  For  t  >  Tmax,  ionization  processes  become  significant  and 
then  the  first  term  on  the  right  hand  side  of  Eq.  (13)  is  no  longer  negligible  in  comparison 
to  the  second  one.  Finally,  Fig.  3  clearly  shows  that  a  quasi-steady-state  condition  has 
to  be  reached  to  determine  Ay.  In  the  case  of  an  initial  ionizing  situation,  the  general 
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characteristics  of  the  results  are  the  same. 

Finally,  to  determine  kr,  respectively  ki,  we  start  from  an  initial  recombining,  respectively 
ionizing  situation,  and,  stop  the  calculations  when  \ne\  and  | iVi  +  A^2  +  A^3 1  agree  to  within 
1%.  Then  we  determine  tqss  and  the  value  of  kr,  respectively  ki,  with  the  instantaneous 
values  of  rie,  and  nl,  respectively  neNa.  The  interest  of  this  approach  is  the  possibility  to 
follow  the  temporal  evolutions  of  populations  and  to  determine  also  Tqss- 
It  is  interesting  to  point  out  that  ionization  and  recombination  rate  coefficients  are  not 
independent.  Indeed,  at  the  stationary  state  of  the  reaction  0+  +  2e_  O  +  e~ ,  the  two 
terms  on  the  right  hand  side  of  Eq.  (13)  have  to  balance  to  give  rie  =  0.  Then,  ki  can  be 
derived  from  kr  by 

kl  =  kr[-0-  =  krKstat  (14) 

We  have  checked  that  both  methods  gave  the  same  results. 


3.4.3  Determination  of  effective  three-body  recombination  rate  coefficients 

Optically  thick  cases 

In  order  to  represent  an  optically  thick  medium,  in  this  section,  all  escape  factors  in 
Eqs.  (8)  and  (9)  are  set  equal  to  zero.  Then,  as  only  elementary  collisional  processes  are 
considered,  the  results  of  this  section  are  denoted  with  the  superscript  C.  It  is  interesting  to 
note  that  these  results  may  be  applied  to  optically  thin  media  where  the  electron  number 
density  is  high  enough  to  neglect  radiative  processes  in  comparison  to  collisions.  After 
having  considered  different  initial  conditions,  we  have  pointed  out  that  in  an  optically 
thick  medium,  the  three-body  recombination  rate  coefficient  for  oxygen  depends  only  on 
the  electron  temperature  (Fig.  5).  We  have  noted  no  influence  of  either  densities  or  initial 
distribution  on  atomic  levels.  The  simplest  best-fit  curve  for  our  computed  result  is 

kcr  =  1.32  x  lir39  f  mV1  (15) 


for  8000  K<  Te  <  20  000  K.  Park  (1990)  assumed  that  the  three-body  recombination  rate 
of  oxygen  is  the  same  as  the  one  of  nitrogen,  derived  for  4000  K<  Te  <  20  000  K  in  using 
for  all  levels  the  elementary  rates  derived  by  Gryzinski  (1965)  for  hydrogen 


krpark  —  1.15  X  10 


(16) 


Figure  5  shows  that,  over  the  whole  temperature  range  of  this  study,  the  value  obtained 
with  our  model  is  on  the  average  five  times  smaller  than  the  one  proposed  by  Park. 
With  our  method,  we  have  also  calculated  the  effective  three-body  recombination  rate 
coefficient  (denoted  k^)  obtained  when  one  uses  for  all  levels  the  elementary  rates  derived 
by  Gryzinski  (1965)  for  hydrogen.  In  this  case,  Fig.  5  shows  that  the  value  of  kf1  is 
close  to  the  one  proposed  by  Park.  Cacciatore  and  Capitclli  (1976)  proposed  to  modify 
the  elementary  rate  coefficients  derived  by  Gryzinski  to  take  into  account  that  low-lying 
levels  of  atomic  oxygen  are  nonhydrogenic.  That  is,  Gryzinski  values  are  multiplied  by 
a  factor  4  for  the  ground  state.  It  is  interesting  to  mention  that  this  correction  has 
no  influence  on  the  calculated  value  of  k l1  for  atomic  oxygen.  Finally,  on  Fig.  5,  the 
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fitting  law  of  the  three-body  recombination  rate  of  nitrogen  obtained  in  a  similar  work 
on  a  collisional-radiative  model  for  atomic  nitrogen  is  represented.  This  law  is  given  for 
4000/i  <  Te  <  12  000/i  and,  therefore  should  be  handle  with  care  for  Te  >  12  000/i. 
However,  it  is  interesting  to  note  that  three-body  recombination  rates  for  oxygen  and 
nitrogen  remain  close  to  each  other  for  <8000 A'  <  Te  <  20  000/d. 


T.[K] 


Figure  5:  Recombination  rate  in  optically 
thick  cases  as  a  function  of  Te.  Solid 
line:  our  model.  Dot-dashed  line:  Park 
(1990).  Dashed  line:  our  atomic  model 
with  Gryzinski  (1965)  cross  sections.  Solid 
line  with  open  circles:  Recombination  rate 
of  nitrogen  (Bourdon  and  Vervisch  (1996)). 


T  [K] 


Figure  6:  Recombination  rate  in  optically 
thick  cases  as  a  function  of  Te.  Thick  solid 
line:  our  40  level  model.  Thin  solid  line: 
9  levels.  Solid  line  with  open  circles:  15 
levels.  Dashed  line:  20  levels.  Dot-dashed 
line:  30  levels. 


Influence  of  the  number  of  levels 

Figure  6  shows  the  influence  of  the  number  of  atomic  levels  on  the  calculated  effective 
three-body  recombination  rate  coefficient  value  in  optically  thick  cases.  As  the  number 
of  atomic  levels  increases,  the  calculated  recombination  rate  converges  towards  an  upper 
limit  and  the  lower  the  electron  temperature  is,  the  more  numerous  atomic  levels  are  re¬ 
quired  to  determine  k^ .  This  figure  puts  forward  clearly  that  the  9-level  model  proposed 
by  Soon  and  Kune  (1990)  is  unadapted  to  our  temperature  range  and  that  it  was  nec¬ 
essary  to  consider  a  more  sophisticated  atomic  model.  Our  forty  level  model  appears  to 
be  a  good  compromise  between  the  inaccuracy  of  atomic  models  for  highly  excited  levels 
and  the  reliability  of  the  effective  rate  coefficients  derived  in  this  study.  It  is  important 
to  mention  that  the  determination  of  at  low  electron  temperatures  Te  <  2  000  K  is 
generally  difficult  since  it  depends  strongly  on  highly  excited  levels.  We  have  noted  in 
this  case,  that  the  more  numerous  highly  excited  atomic  levels  are  taken  into  account,  the 
higher  the  derived  recombination  rate  is. 

Optically  thin  cases 

In  optically  thin  cases,  the  results  obtained  for  different  initial  recombining  conditions  in¬ 
dicate  that  the  three-body  recombination  rate  depends  only  on  the  electron  temperature 
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and  the  electron  number  density.  First,  neglecting  radiative  and  dielectronic  recombina¬ 
tions,  Fig.  7  shows  that  spontaneous  emission  increases  the  value  of  the  recombination 
rate  (denoted,  in  this  case,  kEE).  As  expected,  at  high  electron  number  densities,  the 
recombination  rate  converges  to  the  value  obtained  in  optically  thick  cases. 


TJK] 


Figure  7:  Influence  of  spontaneous  emis¬ 
sion  on  the  recombination  rate  in  an  opti¬ 
cally  thin  plasma.  Numbers  at  right  denote 
x  with  ne  =  10x  m~3.  The  thick  solid  line 
corresponds  to  optically  thick  cases. 


Figure  8:  Recombination  rate  derived  in 
taking  into  account  all  elementary  pro¬ 
cesses,  in  an  optically  thin  plasma.  Same 
notations  as  Fig.  7.  The  thin  solid  line 
corresponds  to  the  expression  proposed  by 
Park  (1990). 


Conversely,  for  low  electron  number  densities,  kEE  tends  towards  an  upper  limit.  If  all 
the  elementary  processes  of  section  3.2  are  considered,  Fig.  8  shows  that  the  influence 
of  radiative  and  dielectronic  recombinations  increases  as  ne  decreases.  In  this  case,  the 
recombination  rate  has  no  upper  limit  at  low  electron  number  densities.  These  results  are 
in  qualitative  agreement  with  those  obtained  earlier  by  Taylor  and  Ah  (1986).  Figure  8 
indicates  that  Park’s  expression  [Eq.  (16)],  initially  proposed  for  an  optically  thick  case, 
corresponds  to  an  optically  thin  plasma  where  1020m-3  <  ne  <  1021m-3  and  8000A'  < 
Te  <  20  000/i . 

In  order  to  implement  the  results  of  this  work  in  numerical  studies  on  air  plasma  flows, 
simple  analytical  laws  are  required.  In  the  general  Eq.  (9),  radiative  and  dielectronic 
terms  can  be  isolated,  then  Eq.  (13)  can  be  written 


kineNa  -  kEEn3e  -  aRn2e 


(17) 


where  k(r'E  takes  into  account  all  the  elementary  collisional  processes  and  spontaneous 
emission  (Figure  7  ),  and 

aR  =  aRR  +  aRD  (18) 

The  total  effective  recombination  rate  kr  is  therefore  kEE  +  aR /ne.  This  formulation  is 
of  great  interest  since  aRR  [Eq.  (7)]  and  aRD  [Nussbaumer  and  Storey  (1983)]  can  be 
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expressed  analytically.  Then,  we  fitted  kfE  for  different  electron  number  densities  to  the 
following  expression 

kfE  =  ATb ,  mV1  (19) 

for  8000/1  <  Te  <  20  000/1 .  The  coefficients  A  and  b  are  given  in  Table  2.  The  discrepancy 
between  the  total  recombination  rate  value  calculated  with  the  analytical  expressions  [Eqs. 
(19)  and  (18)]  and  the  one  determined  directly  and  represented  on  Fig.  8  remains  within 
a  factor  of  2. 


ne  (m  3) 

A  (m6s  4 5K  b ) 

b 

1018 

7.27  x  10"31 

-1.73 

1019 

4.29  x  10"30 

-1.95 

102° 

2.58  x  10"28 

-2.51 

1021 

3.19  x  10-25 

-3.43 

Table  2:  A  and  b  parameters  for  the  effective  three-body  recombination  rate  coefficient 
kfE  =  ATb,  m6s_1  for  different  electron  number  densities 


Comments  on  the  definition  of  kr 

Following  the  method  described  above  to  determine  aR,  one  could  propose  to  derive  the 
effective  three-body  recombination  rate  coefficient  directly  from  Eq.  9 

K  =  YJPd(Te)  =  f(Te)  (20) 

i 

By  definition,  this  coefficient  depends  only  on  Te  and  takes  into  account  only  direct 
recombination  processes  on  the  atomic  levels,  and  neglects  the  other  processes  which  link 
the  atomic  level  populations  together.  Figure  9  shows  that  this  coefficient  k*  is  slightly 
higher  than  kfEax  which  is  the  maximal  value  of  kfE.  This  coefficient  has  no  physical 
interest  to  describe  the  0+  +  2e_  — »  O  +  e~  global  recombination  process.  This  result 
already  mentioned  by  Park  (1969)  is  clearly  illustrated  here. 

3.4.4  Determination  of  the  effective  ionization  rate  coefficient 

Optically  thick  cases 

In  optically  thick  cases,  the  results  obtained  for  different  initial  ionizing  conditions  indicate 

that  the  effective  ionization  rate  coefficient  (denoted  kf)  depends  only  on  the  electron 
temperature.  This  can  be  simply  derived  from  Eq.  (14),  since  in  this  case,  kf  depends 

only  on  Te  and,  the  stationary  state  of  the  reaction  0+  +  e~  O  +  2e_  corresponds  to  a 

Saha  equilibrium,  therefore  Kstat  =  KSaha(Te ).  We  have  fitted  the  equilibrium  constant 
to  the  following  expression 

5 

ln(Kstat)  =  ^jAiZi,  (21) 

i=0 

where  Z  =  Zn(104/Te),  for  8000/i  <  Te  <  20  000/1.  For  a  Saha  equilibrium,  the  coeffi¬ 
cients  Ai  are  given  in  Table  3.  Figure  10  compares  the  value  of  kf  obtained  in  our  study 
with  three  other  values  often  used  in  high  temperature  air  kinetic  schemes,  but  which 
accuracy  is  unknown: 
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•  The  effective  ionization  rate  coefficient  proposed  by  Losev  et  al.  (1994)  which  is 
derived  from  the  ionization  cross  section  from  the  ground  state  measured  by  Smith 
et  al.  (1962). 


_,o  n  /  — 157980 

kiLosev  =  8.64  x  10  18  Te°'68  exp  ( - — - 

•  The  value  proposed  by  Gupta  et  al.  (1990) 

_9qi  /  — 158120  \  o 

kiGupta  =  59.8  Te  2-91exp  i - — - J  m3s 

•  The  value  proposed  by  Park  (1990) 

hpart  =  6478.4  Te  3  7sexp  ( — - - j  m3s  : 


(22) 


(23) 


(24) 


Te[K] 


Figure  9:  Different  definitions  of  the  three- 
body  recombination  rate  coefficient,  k* 
is  the  one-way  coefficient  and  k^Ein  and 
kr.max  are  respectively  the  minimal  and 
maximal  values  of  kEE . 


Te[K] 


Figure  10:  Ionization  rate  coefficient  in  op¬ 
tically  thick  cases.  Solid  line:  our  model. 
Dotted  line:  Park  (1990)  Dot-dashed  line: 
Gupta  et  al.  (1990)  Dashed  line:  Losev 
et  al.  (1994). 


First  of  all,  we  note  that  there  are  strong  discrepancies  between  the  different  models.  As 
expected,  the  value  obtained  with  our  model  is  greater  than  the  one  proposed  by  Losev  et 
al.  since  the  latter  takes  into  account  only  the  ionization  process  from  the  ground  state. 
Figure  10  shows  that  the  value  proposed  by  Gupta  et  al.  is  overestimated  by  more  than 
one  order  of  magnitude.  In  fact,  the  value  of  kf  obtained  with  our  model  turns  out  to 
be  in  close  agreement  with  the  one  proposed  by  Park.  This  result  seems  to  be  in  dis¬ 
agreement  with  the  discrepancy  observed  on  Fig.  5  between  our  value  of  the  three-body 
recombination  rate  in  optically  thick  cases  ( kE )  and  the  one  proposed  by  Park  ( krpark )• 
In  fact,  we  have  noted  that  the  ionization  and  recombination  rate  coefficients  proposed 
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IT 

i 

CO 

Ao 

Ai 

A2 

A3 

a4 

A5 

1017 

30.5304 

-34.0190 

-23.6944 

6.18581 

-0.916160 

-17.5735 

1018 

32.8715 

-33.1659 

-22.8399 

7.45801 

5.02374 

-11.8446 

1019 

35.0493 

-31.4172 

-24.8850 

-0.795804 

17.3478 

11.0059 

102° 

36.6971 

-27.7121 

-6.36738 

-49.7549 

-187.362 

-149.453 

1021 

38.6205 

-35.9083 

6.47223 

62.7877 

-35.6484 

-97.7333 

1022 

41.9019 

-45.7554 

-37.4442 

166.309 

443.193 

298.333 

1023 

45.2287 

-31.9500 

-57.7671 

-69.9729 

-17.6998 

19.8520 

j^Saha 

47.0661 

-17.1846 

-7.99356 

-2.64241 

-0.579718 

0.182226 

Tabic  3:  Parameters  in  KSaha  and  in  the  equilibrium  constant  Kstat  for  different  total 
particle  number  densities. 


by  Park  in  optically  thick  cases  are  not  related  by  the  Saha  equilibrium  constant.  That 
is,  krparkKSaha  is  on  the  average  4.7  times  greater  than  kipark\  which  is  also  the  average 
discrepancy  observed  on  Fig.  5  between  k ,,c  and  krpark- 

Optically  thin  cases 

In  optically  thin  cases,  radiative  processes  enhance  recombination  and  then  the  ionization 
degree  at  the  stationary  state  is  lower  than  the  one  corresponding  to  the  Saha  equilibrium. 
Consequently,  Kstat  <  KSaha(Te),  where  Kstat  depends  on  Te  but  also  on  densities.  As 
explained  by  Kune  and  Soon  (1989)  the  most  significative  density  of  the  nonequilibrium 
properties  of  the  stationary  state  is  the  total  particle  density  Nt  (note  that  the  definition  of 
Nt  proposed  by  these  authors  is  the  double  of  ours)  and  not  the  electron  number  density. 
As  expected,  as  the  density  increases,  collisional  processes  become  more  efficient  than 
radiative  ones  and,  then  the  equilibrium  constant  converges  toward  the  Saha  equilibrium 
constant,  all  the  faster  since  the  electron  temperature  is  high.  First,  neglecting  radiative 
and  dielectronic  recombination  processes,  the  equilibrium  constant  tends  toward  a  lower 
limit  for  low  values  of  Nt.  When  radiative  and  dielectronic  recombination  processes  are 
taken  into  account,  this  lower  limit  disappears.  As  Soon  and  Kune  (1990)  considered 
more  accurately  the  reabsorption  of  the  emitted  radiation,  a  direct  comparison  with  their 
results  at  the  stationary  state  is  impossible.  However,  a  qualitative  agreement  is  observed. 
Finally,  to  implement  our  results  in  flow  codes,  we  have  fitted  Kstat  for  different  values 
of  Nt  to  the  expression  of  Eq.  (21)  for  8000A'  <Te<  20  000/1.  The  coefficients  A;  are 
given  in  Table  3.  The  discrepancy  between  the  ionization  rate  coefficient  value  calculated 
with  the  analytical  expression  [Eq.  (21)]  and  the  one  determined  directly  remains  within 
a  factor  of  2  (Bourdon  et  al.  (1998)). 

3.4.5  Study  of  the  quasi-steady-state  relaxation  time 

In  the  previous  sections,  we  have  put  forward  that  effective  three-body  recombination  and 
ionization  rate  coefficients  exist  only  when  the  system  is  in  a  quasi-steady-state.  There¬ 
fore,  before  implementing  these  rate  coefficients  in  flow  codes,  one  has  to  check  the  validity 
of  the  QSS  hypothesis  [Eq.  (11)].  The  time-dependent  approach  used  in  this  work  allows 
us  to  determine  the  relaxation  time  necessary  to  reach  a  quasi-steady-state,  starting  at 
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t  —  Os  from  a  given  initial  condition. 

Optically  thick  cases 

For  optically  thick  cases,  the  results  obtained  for  different  initial  conditions  (either  ioniz¬ 
ing  or  recombining)  indicate  that  the  quasi-steady-state  relaxation  time  (denoted  TqSS) 
depends  only  on  the  electron  temperature  and  electron  number  density  but  not  on  the 
initial  distribution  on  the  atomic  levels.  Figure  11  shows  that  TqSS  depends  only  weakly 
on  Te  but  increases  as  ne  decreases.  For  ne  =  1019m-3,  the  relaxation  time  derived  in 
using  Gryzinski’s  cross  sections  for  all  levels  is  also  represented.  The  latter  is  one  order 
magnitude  shorter  than  the  one  obtained  with  our  model.  It  is  interesting  to  note  that  re¬ 
laxation  times  obtained  in  oxygen  are  very  close  to  those  calculated  in  nitrogen  (Bourdon 
and  Vervisch  (1996)).  Finally,  Fig.  11  shows  that  the  effective  ionization  and  recom¬ 
bination  rate  coefficients  determined  in  this  work  are  adapted  to  study  optically  thick 
plasma  flows,  where  the  characteristic  flow  time  is  about  10-6  s,  only  if  ne  >  10 20m~3. 
For  lower  electron  number  densities  or  shorter  flow  times,  it  is  necessary  to  consider  one 
conservation  equation  for  each  atomic  level  i  of  the  oxygen  atom. 


TJX] 


Figure  11:  Relaxation  time  in  optically 
thick  cases.  Thin  solid  line  with  open  cir¬ 
cles:  ne  =  1018m~3.  Thin  dashed  line: 
ne  =  1019m-3.  Thin  dot-dashed  line:  ne  = 
102°m~3.  Thin  solid  line:  ne  =  1021m~3. 
The  thick  dashed  line  represents  the  results 
obtained  for  ne  =  1019m~3  in  using  the 
cross  sections  derived  by  Gryzinski  (1965). 


TJK] 


Figure  12:  Relaxation  time  in  optically 
thin  cases  with  at  t  =0s,  bi  =  Aj /N*  =1, 
for  a  given  electron  number  density.  Solid 
line  with  open  circles:  ne  =  1018m-3. 
Dashed  line:  ne  =  1019m~3.  Dot-dashed 
line:  ne  =  102°m-3.  Solid  line:  ne  = 
1021m~3. 


Optically  thin  cases 

In  an  optically  thin  plasma,  the  relaxation  time  depends  on  densities  and  Te  but  also  on 
the  initial  distribution  on  atomic  levels.  Therefore,  it  is  more  difficult  to  derive  general 
results.  However,  we  have  noted  that  for  initial  recombining  conditions  with  at  t  =0  s 
bi  =  Ni/N*  <  1  (where  N*  is  the  Saha  density  of  level  i)  for  a  given  electron  number 
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density,  the  QSS  relaxation  time  depends  only  on  ne  and  Te.  In  these  cases,  Fig.  12 
shows  that  the  QSS  relaxation  time  depends  weakly  on  Te  and  increases  as  ne  decreases. 
In  comparison  to  an  optically  thick  case,  we  note  that  relaxation  times  are  on  the  average 
one  order  of  magnitude  shorter. 

3.5  Conclusions 

In  this  work,  we  have  shown  how  to  use  a  detailed  collisional  radiative  model  to  determine 
the  effective  ionization  and  three-body  recombination  rate  coefficients  of  oxygen  in  a 
plasma  where  1018m-3  <  ne  <  1021m~3  and  8000/i  <  Te  <  20  000/1 .  The  time-dependent 
approach  used  in  this  work  has  clearly  put  forward  the  necessity  that  the  system  be  in  a 
quasi-steady  state  to  derive  meaningful  rate  coefficients.  This  has  led  us  to  study  more 
precisely  the  relaxation  time  tqss  necessary  to  reach  a  quasi-steady-state. 

In  this  work,  the  treatment  of  radiation  has  been  reduced  to  optically  thin  or  thick  cases. 
It  would  be  interesting  to  take  into  account  a  more  accurate  modeling  of  radiation  escape 
factors  in  order  to  study  more  precisely  the  influence  of  the  reabsorption  of  radiation. 
However,  this  effect  will  not  change  the  general  characteristics  of  the  results  obtained  in 
this  paper. 

4  Collisional-radiative  model  for  air 

The  previous  section  has  shown  that  a  temporal  description  of  the  chemistry  of  a  plasma 
is  the  most  relevant  approach  to  understand  its  behavior  by  a  fundamental  point  of  view: 
this  approach  is  the  only  one  allowing  to  check  the  conditions  for  using  a  ionization  or 
recombination  rate  coefficient.  As  a  result,  we  have  chosen  to  reproduce  such  a  study 
but  for  a  more  complicated  medium  where  molecules  play  a  role  in  order  to  see  if  a  quasi 
steady-state  still  exists.  Moreover,  this  study  is  also  motivated  by  our  main  objective: 
the  more  realistic  calculation  of  the  interaction  between  a  plasma  formed  during  an  at¬ 
mospheric  re-entry  of  a  spacecraft  and  its  thermal  shield.  In  that  case,  the  flow  is  indeed 
very  rapid:  as  illustrated  in  the  following,  the  chemistry  is  therefore  time-dependent  in 
Lagrangian  description. 

4.1  Atomic  and  molecular  model 

In  this  work,  air  considered  is  a  mixture  of  heated  nitrogen  and  oxygen  and  their  products: 
argon  as  other  gas  traces  are  neglected.  Due  to  the  wide  range  of  temperatures  and  pres¬ 
sure  studied,  we  take  into  account  N2,  02,  N,  O,  NO,  N2,  02,  N+,  0+ ,  NO+ ,  02, 0~ 
and  electrons.  For  each  species,  we  consider  electronic  excited  states  except  for  A^+,  0+ , 
0~  and  02  as  shown  in  Table  4.  The  energetic  data  are  due  to  Herzberg  (1950),  Capitelli 
et  al.  (2000a)  and  Park  (1990).  46  levels  for  N  and  the  previous  40  levels  for  O  are  taken 
into  account  in  order  to  calculate  correctly  the  recombination  of  N+  and  0+  respectively 
by  three-body  processes  and  in  order  to  estimate  thoroughly  the  net  population  density 
of  the  metastable  states  2D°  and  2P°  for  nitrogen  and  VD  and  1S  for  oxygen  resulting 
from  the  electron  and  heavy  particles  induced  processes.  In  addition,  that  provides  the 
calculation  of  the  excited  population  densities  which  are  the  source  of  the  radiative  signa¬ 
ture  of  the  medium.  The  levels  considered  for  N  are  those  already  mentioned  by  Bourdon 
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Type 

Species 

State 

n2 

X1E+,  A3E+,  B3nq,  c3nu 

Molecules 

o2 

X3E~,  a1  Ag, 

6'EJ.  B3E- 

NO 

X2n,  A2E+ 

,  b2 n,  c2 n,  b’2 a 

N+ 

.Y'2E+,  A2 n„,  B2 E+,  C2E+ 

Molecular 

ot 

x2nq,  a4 

n„,  A2nu,  54s- 

ions 

NO+ 

xl^i=o^uy  A4n 

0~2 

X  U9 

Atoms 

N 

4S0,.. 

.  (46  levels) 

0 

3P,.. 

.  (40  levels) 

Atomic 

N+ 

3p 

ions 

0+ 

Ago 

0~ 

2P 

Table  4:  Species  taken  into  account  in  the  present  CR  model. 


and  Vervisch  (1996). 

Moreover,  since  the  rate  coefficient  for  dissociative  recombination  of  NO+  has  been  calcu¬ 
lated  for  each  excited  vibrational  level  between  v  =  0  and  v  =  14  (Motapon  et  ah  (2006)), 
the  vibrational  excitation  of  this  ion  is  explicitly  considered  contrary  to  the  case  of  the 
other  molecular  species  for  which  the  assumption  of  equilibrium  at  Te  is  done. 

All  in  all,  this  atomic  and  molecular  model  represents  131  different  levels  and  their  cou¬ 
pling  in  the  frame  of  the  present  physical  situation  is  ensured  by  the  following  elementary 
processes. 


4.2  Elementary  processes 

4.2.1  Radiative  processes 

We  take  into  account  45  spontaneous  emission  lines  for  N  and  24  for  O  for  which  Einstein’s 
probability  is  higher  than  5  x  103  s-1  for  N  and  3  x  105  s_1  pour  O.  In  addition,  we  account 
for  the  radiative  and  dielectronic  recombinations.  The  calculations  and  data  for  all  these 
radiative  processes  are  similar  to  those  of  the  previous  section. 

The  possible  re-absorption  of  the  radiations  is  estimated  by  the  following  escape  factor: 


T{1) 


f1  e  k°lxdx 

/o  v/ — 7r  lnx 


(25) 


derived  from  the  work  of  Holstein  (1947)  with  a  typical  plasma  thickness  equal  to  that  of 
the  boundary  layer  (typically  of  the  order  of  l  =  1  cm)  and  assuming  a  Doppler  broadening 
for  which  k o  is  the  absorption  coefficient  by  unit  length  at  the  centerline.  Any  pressure 
broadening  is  assumed  negligible,  the  maximum  pressure  level  being  atmospheric.  In  the 
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same  way,  the  role  of  charged  particles  is  not  considered  since  the  electron  density  ne  is 
sufficiently  weak. 

Some  of  the  states  of  the  species  mentioned  in  Table  4  radiate  importantly.  The  (3  (B2U  — » 
X2n)  and  7  (H2£+  — >  X2n)  systems  of  NO  as  well  as  the  first  positive  (B3 II9  — ■>  H3£+) 
and  the  second  positive  (C3IIU  —*■  B3Ug)  systems  of  N2  and  the  hrst  negative  system 
(_B2£+  — >  X2£+)  of  fV^  have  been  considered.  Since  a  vibrational  equilibrium  is  assumed 
for  these  species  at  the  electron  temperature  Te,  the  equivalent  transition  probability  for 
the  latter  systems  is  calculated  by: 


A 


Eu'(E?;”  Av'v")e 


e" 


G(v') 


G(v’) 

kBTe 


(26) 


where  G{y')  is  the  vibrational  energy,  Av>v«  the  vibrational  transition  probability  and  ks 
the  Boltzmann  constant.  Using  the  data  calculated  by  Laux  and  Kruger  (1992),  we  have 
fitted  each  value  of  A  by: 

3 

A  =  J2a  T  <27) 

i= 0 

whose  coefficients  a,  are  calculated  for  2000  K  <Te<  12000  K  and  may  be  obtained  on 
demand. 


4.2.2  Collisional  processes 

Electron  induced  processes 

When  available,  the  cross  sections  are  chosen  in  accordance  with  the  more  recent  experi¬ 
mental  data.  Since  the  Maxwell  equilibrium  assumption  is  done  for  most  of  the  species, 
the  cross  sections  have  been  integrated  and  the  rate  coefficient  obtained  finally  fitted 
under  the  form: 

kl^{Te)  =  AT;ae-6^  (28) 

In  Table  5  are  summed  up  the  parameters  A,  a  and  6  hence  calculated  for  2000  K  < 
Te  <  10000  K .  The  rate  coefficient  for  the  excitation  of  molecular  species  under  electron 
impact  has  been  calculated  recently  by  Tculet  et  al.  (1999)  under  the  form  of  Eq.(28)  and 
are  used  in  the  present  work. 

There  are  few  data  concerning  the  ionization  of  N  and  O,  especially  for  the  highly  excited 
states.  For  the  ground  state,  we  have  used  for  both  atoms  the  cross  sections  reviewed 
by  Tawara  and  Kato  (1999).  Conversely,  we  have  preferred  the  cross  sections  obtained 
by  Kim  and  Desclaux  (2002)  for  metastable  nitrogen  and  assumed  that  the  cross  section 
for  metastable  oxygen  has  the  same  form  when  the  shift  due  to  the  different  threshold  is 
taken  into  account.  The  derived  rate  coefficients  expressed  by  Eq.(28)  use  the  parameters 
mentioned  in  Table  5.  For  higher  excited  levels,  the  classical  cross  sections  of  Drawin  are 
used  considering  that  each  atom  has  only  one  optical  electron. 

The  data  related  to  the  ionization  of  molecules  are  taken  from  the  work  of  Tculet  et  al. 
(1999)  as  those  concerning  the  dissociation. 

Heavy  particles  induced  processes 

The  excitation  of  atoms  due  to  atom  collisions  has  been  discussed  by  Capitelli  et  al. 
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i 

J 

A 

a 

9 

Ref. 

N(AS°) 

N(2D ) 

1.720  x  RT15 

-0.280 

28250 

Berrington  et  al.  (1975) 

1V(4S'°) 

N(2P ) 

4.004  x  10“16 

-0.324 

42860 

Berrington  et  al.  (1975) 

0{3P ) 

0{lD) 

1.374  x  10”17 

-0.566 

24600 

Itikawa  et  al.  (1990) 

0(3P ) 

OfS) 

9.352  x  10~19 

-0.620 

50110 

Itikawa  et  al.  (1990) 

N^S0) 

N+(3P) 

8.583  x  10”16 

-0.276 

166870 

Tawara  and  Kato  (1999) 

N(2D) 

N+(3P ) 

7.732  x  10~18 

-0.711 

139900 

Kim  and  Desclaux  (2002) 

N(2P ) 

N+(3P ) 

1.026  x  10”17 

-0.677 

126310 

Kim  and  Desclaux  (2002) 

0(3P ) 

0+(4S°) 

5.393  x  10~18 

-0.750 

158330 

Tawara  and  Kato  (1999) 

0{lD) 

0+(4^°) 

8.571  x  10~18 

-0.695 

134190 

Kim  and  Desclaux  (2002) 

O^S) 

0+(4^°) 

1.397  x  10”17 

-0.647 

108800 

Kim  and  Desclaux  (2002) 

Tabic  5:  Parameters  A,  a  and  9  for  the  calculation  of  the  rate  coefficient  (in  m3  s  4)  given 
by  Eq.(28)  for  the  transition  i  —>  j  induced  by  electron  collision.  6  is  expressed  in  K. 


(2000a)  between  the  ground  and  metastable  levels.  We  have  consequently  considered  all 
these  data  in  our  CR  model.  On  the  other  hand,  we  have  adopted  a  linear  form  for  the 
cross  section  crd  with  the  relative  kinetic  energy  e  of  the  colliding  particles  for  excitation 
from  the  ground  or  the  metastable  states  towards  highly  excited  levels: 

4~Aie)  =  ^-A{e-Eii)  (29) 

where  E,Jt  is  the  threshold.  Since  the  parameter  (3Aj~A  weakly  depends  on  the  mass 
of  the  colliding  particles  (Drawin  and  Emard  (1973)),  we  have  assumed  that  the  cross 
section  can  be  derived  from  the  case  of  argon  (Bultel  et  al.  (2002))  taking  the  mean 
value  /3^ =  3.4  x  10~25  m2  eK_1.  For  excitation  between  higher  excited  levels,  we  have 
considered,  as  in  the  case  of  argon,  that  Eq.(29)  is  applicable  with: 

fiA~A  =  8  69  x  10— 22^-2.26  777,2  eyl  (30) 


The  cross  section  for  excitation  of  atoms  under  molecule  impact  behaves  approximately 
as  (Lotz  (1968)): 


cr 


A—M  | 
ij  ' 


e)  =  a0 


In  (c/Ejj) 
e/ Eji 


(31) 


where  ao  —  10~2°m2.  This  form  has  been  adopted. 

Equation  (31)  is  also  assumed  to  condition  the  excitation  of  any  molecular  species  under 
atom  or  molecule  impact  except  in  the  case  when  experimental  data  exist.  We  have  then 
preferred  the  rate  coefficients  compiled  by  Teulct  et  al.  (2001),  Capitelli  et  al.  (2000a) 
and  Kossyi  et  al.  (1992). 

From  these  sets  of  rate  coefficients,  we  have  also  extracted  and  used  those  related  to 
charge  transfer  (with  possible  reassociation  of  atoms  to  form  a  molecule),  reassociation, 
excitation  transfer,  ionization,  dissociation,  neutralization  of  0~  and  O 7  and  attachment 
and  detachment.  Due  to  the  pressure  levels  involved  in  our  CR  model,  we  have  also  taken 
into  account  the  previous  processes  occurring  when  a  third  particle  interacts. 
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Two  reactions  are  known  to  highly  influence  the  distribution  of  nitrogen  and  oxygen 
between  atomic  and  molecular  systems.  These  reactions  called  Zeldovich’s  reactions  are: 

02  +  N  NO  +  O  (32a) 

N2  +  O  Q  NO  +  N  (32b) 

Recently,  Bose  and  Candler  have  up-dated  their  rate  coefficient  using  a  quasiclassical 
trajectory  method  performed  starting  from  an  ab  initio  potential  surfaces  calculation. 
They  have  obtained: 

fc!(TA)  =  4.14  x  10-21  T1/79  e-4000'5/^  m 3  s”1  (33) 

for  the  process  (32a)  (Bose  and  Candler  (1997))  and: 

k2(TA )  =  9.45  x  1(T 18  T9/29  e~A29A9/TA  m 3  s'1  (34) 

for  the  process  (32b)  (Bose  and  Candler  (1996)).  The  rate  coefficients  (33)  and  (34)  are 
used  in  our  work. 

The  processes  (32a)  and  (32b)  can  contribute  to  the  destruction  of  02  and  N2  and  forma¬ 
tion  of  NO.  Of  course,  the  direct  dissociation  and  the  reverse  process  by  heavy  particles 
impact  have  to  be  taken  into  account.  The  rate  coefficients  used  are  due  to  Losev  and 
Shatalov  (1990)  who  have  assumed  the  vibrational  mode  in  equilibrium  with  the  electrons 
at  Te.  We  refer  to  their  paper  for  more  details. 

4.2.3  Dissociative  recombination  of  N%,  and  NO+ 

The  dissociative  recombination  (DR)  of  the  molecular  ions,  which  is  the  subject  of  a  very 
intensive  work  (Guberman  (2003)),  is  well-known  to  play  a  very  important  role  in  the 
case  of  recombining  plasmas  even  at  low  pressure  (Bultel  et  al.  (2002)).  In  addition,  the 
inverse  process,  the  associative  ionization,  allows  in  much  cases  the  formation  of  the  first 
electrons  in  shock  tube  (Drawin  and  Emard  (1973))  as  well  as  for  the  re-entries  problems 
and  consequently  explain  the  ionizing  situations.  In  our  case,  since  N^,  and  NO+ 
are  present  in  the  plasmas  described  here,  dissociative  recombination  has  to  be  considered. 

Dissociative  recombination  of  N. £ 

The  dissociative  recombination  of  has  been  investigated  in  the  CRYRING  storage  ring 
by  Peterson  et  al.  (1998)  over  a  wide  range  for  electron  energy  (10  meV  <  e  <  30  eV). 
They  have  shown  that  its  rate  coefficient  is  weakly  dependent  on  the  vibrational  state  v 
of  the  molecular  ion  and  can  be  set  equal  to: 

aN}  =  1.75  x  10-13  f  ^  J  mV1  (35) 

In  addition,  they  have  obtained  the  branching  fractions  ^(N^)  for  v  =  0  at  OeV  collision 
energy  for  the  following  products  formation: 

N+  +  e~  }  N(4S°)  +  N(2D)  (36a) 
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N+  +  e~  72-^?}  N(4S°)  +  N(2P)  (36b) 

N+  +  e-  73iS}  N(2D)  +  N(2D)  (36c) 

with  7i(A'r21')  =  0.37,  ^2(^2)  =  0.11  and  73 (A^)  =  0.52.  The  collision  energy  range  of 
their  experiment  includes  our  own  (derived  from  the  electron  temperature  Te  levels  con¬ 
sidered  in  this  work).  Moreover,  although  the  rate  coefficient  found  weakly  depends  on 
the  vibrational  excitation  of  the  molecular  ion,  the  vibrational  temperature  of  in  the 
experiment  of  Peterson  et  al.  corresponds  well  to  the  case  of  the  plasmas  considered  in 
this  paper  (Tv  =  Te  ~  6000  K  on  average).  For  all  these  reasons,  their  results  have  been 
integrated  in  our  calculations. 

Dissociative  recombination  of  Ot 

The  dissociative  recombination  of  0%  has  been  investigated  in  the  same  heavy  ion  storage 
ring  by  Peverall  et  al.  (2001)  but  only  on  the  collision  energy  range  1  meV  <  e  <  3eV 
due  to  the  collapse  of  the  cross  section  for  higher  energies.  The  rate  coefficient  derived  is: 

a  +  =  2.4  x  10-13  (  —  )  m3s_1  (37) 

Vsooy  v  ; 

where  the  molecular  ion  is  in  the  vibrational  fundamental  state.  The  branching  fractions 
ascribed  to  the  processes: 


+  0(3P)  (38a) 

+  0(lD )  (38b) 

+  Oi}D)  (38c) 

are  respectively  equal  to  0.20,  0.40  and  0.40.  Although  the  branching  fraction  for  the 
process: 

0+  +  e~  ->  0(lD)  +  OCS )  (39) 

is  not  equal  to  0,  its  order  of  magnitude  and  the  fact  that  it  becomes  almost  0  for  colli¬ 
sion  energies  between  10  and  40  meV  by  increasing  up  to  0.06  for  e  =  300  meV  led  us  to 
consider  finally  that  this  channel  is  closed  and  that  the  dissociative  recombination  for  O 2 
occurs  by  the  processes  (38a),  (38b)  and  (38c)  only.  Founding  no  more  information  in  the 
literature  about  a  possible  influence  of  the  vibrational  excitation  of  O 2  on  the  dissociative 
recombination  rate  coefficient,  we  have  assumed  no  influence  and  consequently  adopted 
Eq.(37)  and  the  previous  branching  fractions. 

Dissociative  recombination  of  NO+ 

The  case  of  NO+  has  to  be  particularly  discussed:  as  we  will  see  in  the  subsequent  part 
of  this  paper,  this  molecular  ion  plays  a  very  important  role.  As  in  the  previous  case 
of  and  0£,  three  factors  condition  the  discussion:  the  vibrational  excitation  of  the 
molecular  ion,  the  collision  energy  range  and  the  branching  fractions.  Early  works  were 
devoted  to  the  experimental  determination  of  the  rate  coefficient  in  room  temperature 
conditions  for  kinetics  in  post- discharge  purposes  (Gunton  and  Shaw  (1965))  (Weller  and 


0+  +  e“  0(3P ) 

0+  +  e-  72-H}  0(3P ) 
0+  +  e-  73iS]  OCD) 
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Biondi  (1968)):  NO+  is  then  on  its  vibrational  ground  state  and  the  collision  energy 
range  is  very  weak.  Other  determinations  of  the  rate  coefficient  were  also  performed  for 
larger  collision  energy  ranges  (Huang  et  al.  (1975)).  With  the  development  of  the  hy¬ 
personic  flights,  experiments  were  done  also  in  shock-tubes  where  the  rate  coefficient  was 
determined  (Dunn  and  Lordi  (1969))  or  calculated  (Hansen  (1968)):  the  vibrational  exci¬ 
tation  of  the  molecular  ion  was  higher  consequently  as  well  as  the  collision  energy  range 
(Te  <  6000  K  in  order  of  magnitude).  Several  theoretical  treatments  have  been  applied 
to  explain  the  experimental  values.  We  mention  the  work  of  Bardsley  (1983)  and  those 
of  Sun  and  Nakamura  (1990)  whose  results,  obtained  with  different  methods,  reveal  that 
the  cross  section  for  v  >  0  can  be  higher  by  one  order  of  magnitude  than  for  v  =  0  when 
e  >  1  eV.  Since  a  calculation  of  the  Maxwell  rate  coefficient  needs  the  value  of  the  cross 
section  up  to  e  =  10  eV  for  Te  =  leV,  the  vibrational  excitation  of  NO+  may  have  an 
influence  on  the  global  dissociative  recombination.  In  order  to  use  an  accurate  value  for 
the  cross  section  on  the  previous  collision  energy  range  and  in  particular  in  order  to  access 
to  the  role  of  the  vibrational  excitation,  we  have  used  the  results  of  the  extensive  calcu¬ 
lations  using  the  Multichannel  Quantum  Defect  Theory  (MQDT)  developed  by  Motapon 
et  al.  (2006)  resulting,  among  other  data,  in  the  cross  section  <Jdr(€,v)  of  the  dissociative 
recombination: 

ATO+(X1S+,u)  +  e.-(e)  a°-^v)  N*  +  O*  (40) 

This  kind  of  reaction  can  occur  in  fact  according  to  two  channels. 


(o)  (6) 


Figure  13:  Schematic  view  of  the  two  chan¬ 
nels  allowing  the  dissociative  recombina¬ 
tion  of  a  molecular  ion  AB  with  electron. 
Left:  the  direct  process,  right:  the  indirect 
process.  The  electron  energy  is  e.  AB**  is  a 
dissociative  potential  curve  of  the  molecule 
AB. 


Figure  14:  Cross  section  <jdr{(-,v)  for  dis¬ 
sociative  recombination  of  NO+  with  elec¬ 
tron  for  different  vibrational  excitation  of 
the  ion  (0  <  v  <  5). 


First,  the  DR  can  be  due  to  a  direct  process:  it  leads  to  the  formation  of  a  doubly  excited 
state  of  NO,  denoted  NO**,  which  dissociate  spontaneously  afterwards  by  producing 
separated  N  and  O  atoms  possibly  electronically  excited  (see  Fig. 13).  The  estimation  of 
the  cross  section  needs  consequently  an  accurate  knowledge  of  the  potential  curves  for 
these  states.  The  higher  the  energy  for  the  electron,  the  more  numerous  the  potential 
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curves.  Second,  the  DR  can  also  occur  through  an  indirect  process:  this  process  leads  to 
the  formation  of  a  Rydberg  level  of  NO  associated  with  a  vibrationnally  excited  state  of 
NO+  which  is  predissociated  by  the  NO**  state.  In  this  case,  the  Rydberg  states  being 
numerous  by  definition,  the  calculation  of  the  contribution  of  the  indirect  channel  to  the 
cross  section  needs  a  great  amount  of  molecular  data. 

Both  channels  need  finally  numerous  data  and  one  of  the  main  issues  to  be  addressed  is  the 
thorough  determination  of  the  potential  curves  related  to  the  states  involved.  After  such 
a  determination,  several  methods  can  be  applied  in  order  to  calculate  the  cross  section. 
One  of  the  most  efficient  one  at  the  present  time  is  based  on  the  Multi-Channel  Quantum 
Defect  Theory.  Its  approach  is  based  on  a  description  of  molecular  states  in  which  only 
part  of  the  electronic  hamiltonian  is  diagonalized,  within  subspaces  of  electronic  states 
with  similar  nature.  We  refer  to  the  paper  of  Motapon  et  al.  (2006)  for  more  details. 

We  can  observe  on  Fig.  14  the  behavior  of  the  cross  section  resulting  from  the  energy  of 
the  colliding  electron  and  the  vibrational  excitation  of  the  molecular  ion  calculated  by 
the  MQDT  approach.  The  influence  of  the  great  number  of  states  is  obvious:  the  cross 
section  presents  much  oscillations  which  correspond  to  resonances  when  the  energy  of 
the  electron  is  accorded  with  the  difference  between  two  particular  states  involved.  In 
addition,  the  global  behavior  does  not  follow  a  e_1  law  contrary  to  that  usually  assumed. 
In  particular,  when  the  energy  of  electron  exceeds  7  eV,  the  cross  section  collapse. 

The  Maxwellian  rate  coefficient  of  the  DR  is  derived  from  the  calculation  of  the  cross 
section.  Each  rate  coefficient  is  finally  fitted  under  the  general  modified  Arrhenius  law 
(28).  We  can  find  in  Table  6  the  parameters  A,  a  and  9  for  each  value  of  the  vibrational 
quantum  number  v  between  0  and  5.  The  usual  form  adopted  for  the  cross  section  leads 


V 

A 

a 

9 

0 

8.486  x  10“n 

0.793 

859.99 

1 

3.958  x  10"14 

0.076 

-661.83 

2 

1.664  x  KT11 

0.704 

-39.13 

3 

6.335  x  1(T13 

0.383 

-282.98 

4 

2.970  x  10-13 

0.311 

-632.21 

5 

2.666  x  10-13 

0.297 

-472.45 

Table  6:  Parameters  A,  a  and  9  for  the  calculation  of  the  rate  coefficient  (in  m3  s-1)  given 
by  Eq.(28)  for  the  dissociative  recombination  of  NO+(X1'E,v)  with  electrons  used  in  this 
work  Motapon  et  al.  (2006).  6  is  expressed  in  K . 

to  a  rate  coefficient  following  a  T~0'5  law.  The  characteristic  temperature  9  being  rather 
small,  the  temperature  dependence  is  mainly  due  to  a.  Table  6  shows  that  the  rate 
coefficient  is  relatively  different  from  a  value  of  the  vibrational  quantum  number  v  to 
another  with  a  form  far  from  a  T-0,5  law:  as  a  result,  the  usual  temperature  dependence 
is  consequently  not  realistic  and  too  simplistic.  Assuming  an  equilibrium  between  the 
vibration  of  NO+  and  electrons,  these  values  lead  to  that  of  the  global  rate  coefficient  for 
the  dissociative  recombination: 

kDR(T)  =  1.658  x  io-10T-a933  e“377-68/T  m3s~1  (41) 
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where  the  main  temperature  dependence  is  due  to  the  T~°'933  term.  This  form  reveals 
that  the  common  law  adopted  for  the  DR  is  erroneous.  This  conclusion  can  be  made 
also  in  the  case  of  A^~  and  0\  whose  rate  coefficient  presents  a  temperature  dependence 
according  to  a  T”°"3  and  T~°'  law  respectively  as  experimentally  determined. 

The  more  realistic  view  of  the  process  (40)  needs  to  know  the  excitation  of  nitrogen  and 
oxygen  atoms  after  the  recombination.  As  far  as  we  know,  only  three  works  were  devoted 
to  this  difficult  issue.  The  earliest  one,  due  to  Kley  et  al.  (1977)  and  performed  by 
photoionization  of  buffered  nitric  oxide  by  a  flash  lamp,  reported  that  only  two  channels 
are  open: 

NO+(X1T,+)  +  e~  -^4  N(4S°)  +  0(3P )  (42) 

and 

1V0+(A"1E+)  +  e~  N(2D )  +  0(3P)  (43) 

with  7!  =  0.24  ±  0.06  and  y2  =  0.76  ±  0.06.  It  is  important  to  note  that  these  results 
have  been  obtained  with  a  Maxwellian  distribution  for  the  kinetic  energy  e  of  colliding 
particles.  Controlling  e,  Vej by- Christensen  et  al.  (1998)  have  more  recently  determined 
the  branching  fractions  in  the  ASTRID  ion  storage  ring  for  a  collision  energy  e  =  0,  0.75 
and  1.35  eV.  They  confirm  in  part  the  results  of  Kley  et  al.  but  with  yx  rather  equal 
to  0.15  ±  0.06  for  e  =  OeV  and  not  well  identified  for  e  =  0.75  and  1.35  eV.  For  these 
two  collision  energies,  they  concluded  only  that  the  process  (43)  represents  at  least  65% 
of  the  dissociative  recombination  events.  Finally,  Hellberg  et  al.  (2003)  have  investigated 
the  dissociative  recombination  of  the  X4T,+  and  a3E+  states  in  the  CRYRING  storage 
ring  and  succeeded  for  the  ground  electronic  state  to  measure  the  branching  fractions 
for  e  =  0,1.25  and  5.6  eV  .  Increasing  the  collision  energy  leads  to  a  multiplication  of 
the  branching.  Table  7  illustrates  their  results  compared  to  those  of  Kley  et  al.  and 
Vejby-Christensen  et  al.  Since  our  CR  model  is  elaborated  to  work  over  a  wide  range  of 
temperature,  each  one  corresponding  to  a  Maxwell  distribution  of  the  collision  energy,  the 
temperature  dependent  branching  fraction  7 i(Te)  for  each  i  channel  has  to  be  determined. 
The  values  finally  retained  in  our  work  result  from  the  following  considerations  and  are 
based  on  the  data  of  Hellberg  et  al.  due  to  their  best  precision.  When  the  collision  energy 
is  insufficient,  the  branching  fraction  is  put  equal  0.  This  is  the  case  at  different  energies 
of  all  the  channels  involved  except  those  producing  N(4S°)  +  0(3P)  and  N(2D)  +  0(3P) 
for  which  the  processes  are  exothermic.  Since  the  branching  fraction  for  the  production 
of  Y(4S'°)  +  0(1P)  is  always  0  considering  the  uncertainty  range,  we  have  assumed  that 
this  channel  is  systematically  closed. 

Moreover,  the  branching  fraction  is  assumed  independent  of  the  collision  energy  when 
the  collision  energy  is  greater  than  the  threshold  (which  is  equal  to  0.81,  1.59,  2.78,  3.81 
and  5.0  eV  for  the  processes  3  up  to  7  respectively).  For  the  channels  (42)  and  (43),  the 
values  adopted  are  chosen  over  the  uncertainty  range  in  order  that  the  sum  of  branching 
fractions  be  equal  to  unity.  As  a  result,  the  branching  fractions  obtained  according  to  the 
energy  ranges  are  summed  up  in  Table  8.  The  temperature  dependent  branching  fraction 
7i(Te)  for  the  i-process  is  finally  obtained  by: 


7 i(Te) 


-+oo 


7 i(e)xe  xaDR(x)dx 


xe  x  odr(x)  dx 


(44) 
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ORGANIZATION 


Ref. 

Process 

Process 

Branching 

i 

fraction  7 j 

K.  NO+{y  =  1.5)  +  e"(e  =  0.34  eV)  — >  N(4S°)  ±  0{3P ) 

K.  NO+(v  =  1.5)  +  e_(e  =  0.34  eV)  — >  N(2D )  +  0(3P) 

V.-C.  NO+( 0)  +  e~(0eV)  — *  N{4S°)  ±  0(3P ) 

V.-C.  NO+( 0)  +  e-(Oey)  — *  N(2D )  +  0(3P) 

V.-C.  NO+( 0)  +  e-(0.75,1.35eV)  — >  N(2D)  ±  0(3P) 

H.  NO+(  0)  +  e-(OeV)  — >  N(4S°)  ±  0(3P) 

H.  NO+( 0)  +  e-(0eV)  — *  N(2D)  ±  0(3P) 

H.  NO+( 0)  +  e_(1.25 eV)  — >  7V(4,S0)  +  0(3P) 

H.  NO+l 0)  +  e_(1.25 eV)  — >  N(2D)  +  0(3P) 

H.  VO+(0)  +  e-(1.25eV)  — >  N(4S°)  +  0(1P) 

H.  VO+(0)  +  e-(1.25eV)  — >  N(2P )  +  0(3P) 

H.  VO+(0)  +  e"(5.6eV)  — ♦  N{4S°)  +  0(3P) 

H.  NO+lo)  +  e-(5.6eV)  — *  N(2D)  +  0(3P) 

H.  1VO+(0)  +  e-(5.6eV)  — >  N(2P )  +  0(3P) 

H.  NO+l 0)  +  e-(5.6eV)  — >  N(2D )  +  O^D) 

H.  NO+(0)  +  e-(5.6eV)  — >  N(2P )  +  0(1P) 

H.  NO+(0)  +  e-(5.6eV)  — >  V(2P)  +  O(^) 

H.  NO+(0)  +  e-(5.6eV)  — >  N(2P )  +  0(15) 


1 

2 

1 

2 

2 

1 

2 

1 

2 

3 

1 

2 

3 

4 

5 

6 
7 


0.24  ±0.06 
0.76  ±0.06 
0.15  ±0.06 
0.85  ±0.06 
>  0.65 
0.05  ±0.02 
0.95  ±0.02 
0.10  ±0.06 
0.70  ±0.20 
0.10  ±0.10 
0.10  ±0.07 
0.03  ±0.02 
0.15  ±0.05 
0.11  ±0.05 
0.31  ±0.05 
0.21  ±0.05 
0.09  ±0.05 
0.10  ±0.05 


Table  7:  Comparison  between  experimental  information  available  on  the  branching  frac¬ 
tions  of  the  dissociative  recombination  of  1V0+(X1£+,  v).  The  value  v  =  1.5  means  that 
NO+  is  weakly  vibrationally  excited.  The  value  e  =  0.34  eV  means  that  the  collision 
energy  is  distributed  according  to  a  Maxwell  distribution.  Ah,  V.  —  C.  and  P.  mean  Kley 
et  al.  (1977),  Vejby-Christensen  et  al.  (1998)  and  Hellberg  et  al.  (2003)  respectively. 


Energy 
range  (eV) 

7i 

72 

73 

74 

75 

76 

77 

[0;  0.81[ 

0.07 

0.93 

0.00 

0.00 

0.00 

0.00 

0.00 

[0.81;  1.59[ 

0.15 

0.75 

0.10 

0.00 

0.00 

0.00 

0.00 

[1.59;  2.78[ 

0.03 

0.55 

0.11 

0.31 

0.00 

0.00 

0.00 

[2.78;  3.81  [ 

0.03 

0.34 

0.11 

0.31 

0.21 

0.00 

0.00 

[3.81;  5[ 

0.03 

0.25 

0.11 

0.31 

0.21 

0.09 

0.00 

[5;  10[ 

0.03 

0.15 

0.11 

0.31 

0.21 

0.09 

0.10 

Table  8:  Branching  fractions  7i(e)  for  Aprocess  (see  Table  7)  adopted  for  the  calculation 
of  the  temperature  dependent  branching  fractions  of  the  dissociative  recombination  of 
VO+(X1S+,0)  with  electrons. 
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where  x  =  e/(kBTe)  and  a  on  the  dissociative  recombination  cross  section.  Figure  15 
shows  the  branching  fractions  calculated  with  the  help  of  the  cross  section  determined 
by  the  MQDT  approach.  In  the  following,  we  have  assumed  that  the  branching  fractions 
behave  for  v  >  0  like  for  v  =  0. 


T,(K) 


Figure  15:  Temperature  dependent 

branching  fractions  obtained  with  our  cal¬ 
culated  cross  section  for  the  dissociative  re¬ 
combination  of  NO+(X1T,+,  0)  for  the  dif¬ 
ferent  channels  investigated  by  Hellberg  et 
al.  (see  Table  7  and  8). 


Figure  16:  Considered  situations  of  re¬ 
entry  displaying  the  pressure,  velocity  and 
temperatures  profiles  along  the  stagnation 
streamline. 


The  fact  that  the  rate  coefficient  depends  strongly  upon  v  has  led  us  to  develop  a  state 
to  state  kinetics  for  NO+. 

4.2.4  Vibrational  processes 

Due  to  the  lack  of  experimental  or  theoretical  investigations  allowing  the  determination 
of  the  cross  section  of  vibration- vibration  (VV)  and  vibration-translation  (VT)  processes 
involving  NO+ ,  we  assume  that  this  ion  behave  like  IV2  for  which  Armenise  et  al.  (1996) 
and  Billing  and  Fisher  (1979)  have  produced  very  useful  data.  Conversely,  the  case  of 
vibration-electron  (VE)  processes  like: 

VO+(V1E+,ui)  +  e~  -»■  NO+(X1'E+, vf)  +  e~  (45) 

corresponding  to  inelastic  collisions  when  Vf  >  i\  and  to  super-elastic  collisions  in  the 
reverse  case  can  be  treated  by  the  MQDT  approach.  The  cross  sections  have  been  de¬ 
termined  and  the  Maxwellian  rate  coefficients  derived.  These  rate  coefficients  have  been 
fitted  under  the  form  (28). 
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4.3  Results 

4.3.1  Ionization  situation 


Although  the  CR  model  elaborated  with  the  help  of  the  previous  data  can  work  up  to 
the  atmospheric  pressure,  we  will  discuss  in  the  following  some  typical  re-entry  situations 
only  displayed  for  instance  by  Fig.  16  where  the  pressure  is  clearly  lower. 

We  have  chosen  to  develop  our  study  from  a  temporal  point  of  view  as  already  developed 
in  the  section  devoted  to  the  elaboration  of  CR  model  for  atoms.  This  choice  will  be 
justified  later. 

The  general  balance  equation  for  a  species  X  on  one  of  its  excited  levels  i  can  be  written 
indeed  by  displaying  the  hydrodynamic  derivative  D / Dt  relevant  under  a  Lagrangian 
description  of  the  flow: 


D\X M 

Dt 


+  [X(*)]Vr  =  -V4  + 


\9[xm 


dt 


CR 


(46) 


This  form  is  particularly  well  adapted  to  estimate  the  importance  of  each  contribution 
in  the  balance:  comparing  their  characteristic  time  scale  is  enough.  In  order  to  esti¬ 
mate  correctly  the  one  related  to  the  chemical  source  [d[X(i)\/dt]CR,  we  assume  that  the 
characteristic  time  scale  for  the  fluid  expansion  of  contraction  taken  into  account  in  the 
balance  equation  by  the  term  [X(i)J  V?)  and  that  of  the  diffusion  term  V  Jxl  tend  to 
infinity.  The  calculations  are  performed  assuming  that  the  Maxwell-Boltzmann  equilib¬ 
rium  is  reached  for  the  kinetic  energy  of  atoms  and  molecules  at  Ta  and  electrons  at  Te. 
These  assumptions  are  verified  a  posteriori  by  comparing  the  characteristic  time  scale  for 
reaching  a  Maxwellian  equilibrium  with  that  of  the  chemistry.  Since  the  electrons  are 
produced  owing  to  the  collisions  resulting  from  the  high  temperature  for  heavy  species, 
we  can  have  Te  <T R 

We  can  observe  on  Fig.17  and  Fig. 18  the  time  evolution  of  the  number  density  of  the 
species  of  Table  4  when  a  fluid  particle  with  T  ~  200  K  crosses  the  shock  front  and 
is  suddenly  in  high  temperature  (7000/1  for  Fig.17  and  10000  K  for  Fig. 18)  and  high 
pressure  (6  kPa  for  Fig.17  and  15  kPa  for  Fig.  18)  conditions.  We  can  observe  that  the 
evolution  of  population  densities  starts  immediately  from  t  —  0  and  stops  at  rc  =  3  x  10-5  s 
for  Te  =  Ta  =  7000  K  and  at  rc  =  4  x  10_4s  for  Te  =  5000  K  in  the  case  of  Fig.17.  The 
fact  that  the  evolution  starts  immediately  from  t  =  0  is  due  to  the  very  strong  initial 
non  equilibrium.  The  initial  concentrations  correspond  as  already  mentioned  to  those 
before  the  shock  front  at  high  altitude  typically  where  the  temperature  is  approximately 
T  ~  200  K:  there  is  no  populated  excited  states  in  these  conditions.  As  soon  as  the 
thermodynamic  conditions  change  due  to  the  crossing  of  the  shock  front,  the  density  of 
species  begin  to  evolve.  When  electrons  and  heavy  species  are  in  the  same  thermal  con¬ 
ditions  downstream  the  shock  front,  (Te  =  Ta),  tc  is  shorter  and  the  state  following  this 
moment  is  quite  close  to  chemical  equilibrium,  the  pressure  being  sufficiently  high.  When 
Te  <T a,  this  state  is  also  quite  close  to  chemical  equilibrium  and  is  almost  the  same  than 
the  latter.  One  of  the  important  difference  between  both  cases  we  can  observe  concerns 
tc.  Although  the  electrons  have  a  relatively  weak  influence  on  the  final  state,  their  den¬ 
sity  remaining  low,  decreasing  their  temperature  Te  from  7000  K  to  5000  K  modifies  rc 
by  one  order  of  magnitude.  Figure  17  shows  clearly  that  the  negative  charged  particles 
are  not  electrons  during  the  first  part  of  the  evolution  (t  <  10-6  s  to  10”5  s  according  to 
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the  level  of  Te).  Changing  Te  cannot  therefore  be  ascribed  to  the  action  of  the  electrons. 
This  behavior  is  in  fact  the  result  of  the  dissociation  of  N2,  02  and  NO  under  their  own 
impact  which  leads  to  higher  concentrations  observed  before  10' ' s  for  Te  =  Ta ■  Indeed, 
the  rate  coefficients  that  we  use  are  derived  from  the  paper  of  Losev  and  Shatalov  (1990) 
and  depends  on  both  Ta  and  Te  assuming  an  equilibrium  between  the  electrons  and  the 
vibration  of  the  molecules. 


Figure  17:  Time  evolution  in  a  typical  re¬ 
entry  situation:  ionization  (thermal  equi¬ 
librium:  Te  =  Ta  =  7000  K  in  thick  lines 
and  5000  K  =  Te  <  Ta  =  7000  K  in  thin 
lines  for  a  case  of  thermal  disequilibrium, 
negative  ions  in  dot-dashed  lines). 


N 


Figure  18:  Time  evolution  in  a  typical  re¬ 
entry  situation:  ionization  (thermal  equi¬ 
librium:  Te  =  Ta  =  10000  K  in  thick  lines 
and  8000  K  =  Te  <  Ta  =  10000  77  in  thin 
lines  for  a  case  of  thermal  disequilibrium, 
negative  ions  in  dot-dashed  lines). 


This  hypothesis  is  highly  questionable  in  the  present  case,  the  number  density  of  electrons 
being  very  weak  before  10'7  s:  this  is  confirmed  by  the  evolution  of  Tv°  ^  displayed 
on  Fig.  19.  This  temperature  is  defined  by: 


JlNO+(X) 


1 


kB 


dln{[NO+{X,v)}/ZR{v)) 

dE(X,v) 


Isl 


(47) 


where  ZR{v)  is  the  rotational  partition  function  of  the  vibrational  level  v  considered  whose 
energy  is  E(X,  v )  and  Isl  means  that  the  derivative  is  in  fact  the  slope  of  the  least  squares 

line. 

Before  10'7  s,  the  population  density  of  the  vibrational  states  are  increasing  and  the 
efficiency  of  the  elastic  collisions  with  electrons  or  heavy  particles  is  insufficient  to  ensure 
rapidly  an  equilibrium:  the  vibrational  distribution  does  not  follow  a  Boltzmann  equilib¬ 
rium  neither  with  Te  nor  Ta  and  depends  on  chemical  processes.  This  is  confirmed  by 
Table  9  listing  the  main  processes  responsible  for  the  creation  and  destruction  for  N,  O, 
NO,  O2 , 0~ ,  electrons  and  NO+  for  t  ~  10'7s.  This  list  shows  clearly  that  populating 
the  state  v  =  0  for  A^O+(A"1S+)  is  dominated  by  associative  ionization  between  O  and 
N(2D)  atoms  and  the  inverse  process  of  neutralization  between  NO+  and  O7  in  a  reduced 
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Figure  19:  Time  evolution  of  the  vibrational  temperature  of  NO+(X1H+)  in  the  case  of  Fig. 17, 
the  gas  having  crossed  the  shock  front  and  whose  temperatures  are  such  that  Te  =  Ta  =  7000  K 
(thick  line)  or  5000 K  =  Te  <  Ta  =  7000 IF  (thin  line).  The  thin  dot-dashed  line  refers  to  the 
typical  recombination  inside  the  boundary  layer  with  p  =  2  kPa  from  Te  =  Ta  =  7000  K  to 
Te  =  Ta  =  3000  K  (see  Fig. 20  and  21)  and  the  thick  dot-dashed  line  refers  to  the  same  situation 
but  with  p  =  10  kPa  (see  Fig. 22  and  23). 


Processes 


(«) 


(P) 


->  N 
N  -> 
-►  O 

o  -> 

->  NO 
NO  -»• 


o~ 


0~  ■ 

— >  e 


e  — > 

->  NO+ 
NO+  -> 


N2  +  O  -*•  NO  +  N 
N  +  02  ->  NO  +  OQ-D) 

o2  +  n2^o  +  o  +  n2 

0  +  N2  —  0(1D)  +  1V2 

A^2  +  O  — >  ATO  +  A" 

ATO  +  N2  ^  N  +  O  +  N2 
02(a)  +  O  — ■>  02  +  O 
02  +  02  — ■>  02  +  02  +  e 
02  +  e”  — >  0~  +  O 
O  +  02(a)  — *  02  +  O 
02  +  02  — *  e  +  02  +  02 
e~  +  02  — >  0“  +  O 
O  +  N(2D)  ->  NO+  +  e- 
NO+  +  1V2  -f  7VO+(u  =  1)  +  JV2 


1V2  +  W2  AT  +  AT  +  N2 
N  +  02^  NO +  0 

o2  +  o2  — >0  +  0  +  02 

O  +  N2^  NO  +  N 
02  +  N  ^  NO  +  O^D) 
NO+O^O+O+N 
o2  +  o2  — *  02  +  O2 
02  +0  — >  O  +  02(a) 

02  +0  — >  O  +  02(a) 

O-  +  O  — >  02  +  e“ 

O”  +  0  — >  e_  +  02 
e"  +  0+  -f  O  +  0(1O) 

ATO  +  02  — >  NO+  +  Oi 
NO+  +  02  ->  A^O+(u  =  1)  +  02 


Table  9:  Main  processes  at  time  t  =  10-7  s  during  the  formation  of  the  plasma  illustrated 
by  Fig.  17  when  Te  =  Ta-  The  arrow  indicates  that  the  species  X  is  produced  (— >  X )  or 
destroyed  (X  — *•).  The  second  column  (a)  corresponds  to  the  more  efficient  process  and 
the  third  one  (/?)  to  that  just  follows.  The  species  are  on  their  ground  state  if  there  is  no 
indication  and  NO+  means  NO+(X1T,+,v  =  0)  except  in  the  case  v  —  1. 
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way  while  populating  the  higher  excited  vibrational  states  are  resulting  from  collisions 
with  N2  and  02  molecules.  Afterwards,  elastic  collisions  with  electrons  remain  not  very 
efficient,  ne  remaining  relatively  low:  this  is  illustrated  by  Fig.  17  where  changing  the 
electron  temperature  from  7000  K  to  5000  K  does  not  change  Tv°  ^  by  a  same  factor. 
Anyway,  this  behavior  illustrates  the  high  sensitivity  of  the  kinetics  towards  the  hypoth¬ 
esis  adopted  under  high  degree  of  non  equilibrium  and  the  limitations  of  our  CR  model. 
Table  9  is  also  interesting  because  it  provides  understanding  of  the  global  mechanism  of 
the  plasma  formation  by  coupling  with  Fig.  17.  Although  the  number  density  of  O 2  and 
0~  are  completely  negligible  when  the  steady  state  is  reached,  these  species  play  a  key 
role  in  the  ionization  phase.  It  is  often  assumed  that  the  ionization  in  shock  waves  is 
driven  by  associative  ionization.  When  the  gas  contains  molecules,  they  are  dissociated 
under  heavy  particles  impact  and  the  atoms  produced  associate  to  form  electrons.  This 
is  particularly  true  in  the  case  of  noble  gases  where  atoms  already  exist:  the  associative 
ionization  can  start  immediately  (Drawin  and  Emard  (1973)).  We  see  that  this  is  not  the 
case  for  molecules  being  electronegative  or  having  the  possibility  to  form  such  species. 
The  most  important  negative  charged  particles  are  02  and  0~  whose  formation  is  de¬ 
rived  from  the  charge  transfer  process  02(a)  +  0~  — >  02  +  O  and  the  inverse  process 
of  neutralization  02  +  02  ^  02  +  02  for  02  and  02  +  e~  —>  0~  +  O  and  the  charge 
transfer  process  O.7  +  O  — >  0~  +  02(a)  for  0~.  These  processes  delay  the  formation 
of  free  electrons  and  form  a  reservoir  which  will  be  a  subsequent  source  for  these  parti¬ 
cles:  in  order  of  decreasing  importance,  the  electrons  are  produced  first  by  the  process 
O  +  0~  — >  02  +  e~  and  second  by  the  associative  ionization  O  +  N(2D )  — >  NO+  +  e~  at 
time  t  —  2  x  10-6  s  for  Te  =  T 4.  The  numerical  study  of  air  shock  waves  do  not  usually 
involve  these  processes  (Abe  et  al.  (2002)).  Nevertheless,  it  could  be  interesting  to  adopt 
them  and  observe  their  influence. 

Apart  from  the  problem  of  the  vibrational  distribution  not  necessarily  in  equilibrium  dur¬ 
ing  the  first  part  of  the  evolution,  the  concentrations  calculated  for  t  <  5  x  10-8  s  are 
not  significant  since  it  corresponds  to  too  short  time  scales  in  regard  to  that  for  reaching 
Maxwellian  equilibrium  for  electrons  and  heavy  species.  Moreover,  assuming  a  linear  de¬ 
crease  of  the  velocity  through  the  shock  front  between  and  vA  (see  Fig. 16)  allows  to 
estimate  the  duration  of  crossing  the  shock  front  (having  a  thickness  equal  to  the  mean 
free  path  A  under  upstream  conditions)  of  a  fluid  particle  such  that: 

Tsc  =  - — - In  —  (48) 

Uoo  -  v A  VA 

For  a  mean  free  path  A  ~  1CT3  m  in  order  of  magnitude,  we  have  rsc  ~  10~6  s.  Before  this 
value,  the  fluid  particle  is  not  yet  in  the  conditions  downstream  the  shock  front  and  the 
calculation  has  to  be  considered  prudently.  Nevertheless,  this  calculation  is  interesting 
because  it  provides  the  chemical  scheme  leading  to  the  formation  of  the  plasma  starting 
from  a  cold  gas  when  the  conditions  change  suddenly.  Anyway,  the  state  following  rc  is 
quite  close  to  equilibrium:  it  depends  very  weakly  on  the  initial  conditions.  Assuming 
once  again  a  linear  decrease  of  the  speed  from  vA  to  0  between  y  =  A  and  the  wall,  and 
assuming  also  that  the  cold  gas  is  suddenly  in  the  conditions  just  behind  the  shock  front 
at  t  =  0,  the  y-location  of  the  fluid  particle  at  the  time  t  is: 

y(t)  =  Ae~v*t/A  (49) 
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When  the  near-equilibrium  is  reached  (either  at  the  time  rc  =  3  x  10-5  s  with  Te  =  Ta  = 
7000/1  or  at  the  time  rc  =  4  x  10~4s  for  the  other  case  related  to  Fig.  17),  the  fluid 
particle  is  located  either  in  y  ~  6.3  cm  or  in  y  ~  0.8  cm  respectively.  In  both  cases,  one 
can  conclude  that  the  edge  of  the  boundary  layer  is  therefore  close  to  chemical  equilibrium 
because  y  >  5. 

The  global  behavior  when  Ta  =  10000  K  is  the  same  as  we  can  see  on  Fig. 18.  Due  to 
strong  non  equilibrium,  the  time  evolution  starts  immediately  but  stops  approximately 
at  the  same  time:  when  t  =  rc  =  9  x  10-7  s  for  Te  =  Ta  =  10000  K  and  when  t  = 
rc  =  1.5  x  10-6s  for  8000/1  —  Te  <  Ta  —  10000  K.  This  is  due  to  the  same  processes 
of  dissociation  discussed  previously  whose  rate  coefficient  varies  more  slowly  between 
8000  K  and  10000  K  than  between  6000/1  and  8000/1.  As  in  that  case,  Fig. 18  shows 
that  the  electrons  are  not  the  most  significant  negative  charged  particles.  The  first  part 
of  the  evolution  is  also  questionable  regarding  the  characteristic  time  scales  for  reaching 
Maxwell-Boltzmann  equilibrium  for  electrons  and  heavy  species  as  well  as  regarding  rsc , 
Eq.(48)  leading  to  the  estimation  rac  =  2  x  10”7s.  The  same  conclusion  can  be  derived 
from  Eq.(49):  whatever  the  conditions,  the  fluid  particle  is  not  still  in  the  boundary  layer 
when  the  near-equilibrium  is  reached:  its  location  is  y  ~  4.8  cm  for  thermal  equilibrium 
(Te  =  Ta  =  10000/1)  and  is  y  ~  4.7  cm  when  Te  =  8000/1.  The  chemical  relaxation 
is  therefore  extremely  rapid  and  the  boundary  layer’s  edge  is  therefore  close  to  chemical 
equilibrium.  Finally,  the  most  important  processes  leading  to  the  plasma  formation  are 
globally  the  same  as  those  listed  in  Table  9. 

4.3.2  Recombination  situation 

In  the  literature,  different  simplified  kinetic  schemes  have  been  implemented  in  Navier- 
Stokes  codes  to  study  high  temperature  air  flows  and  particularly  boundary  layers.  Most 
of  these  schemes  were  developed  during  the  70’s  and  80’s  in  the  frame  of  the  shuttle’s 
building.  We  have  chosen  to  discuss  three  widely  used  simplified  kinetic  schemes  due  to 
Dunn  and  Kang  (1973),  Gupta  et  al.  (1990)  and  Park  (1993),  the  latter  being  an  up-dated 
version  of  an  older  one  Park  (1985).  Except  O7  and  0~,  they  consider  the  same  species 
as  those  we  have  taken  into  account.  In  addition,  apart  from  the  scheme  of  Park,  they 
assume  a  thermal  equilibrium  between  electrons  and  heavy  particles. 

Conversely  to  the  case  of  the  kinetic  scheme  of  Park,  those  of  Dunn  &  Kang  and  Gupta 
contain  the  rate  coefficient  of  each  inverse  process:  as  a  result,  we  have  calculated  each  in¬ 
verse  rate  coefficient  for  the  kinetic  scheme  of  Park  starting  from  the  equilibrium  constant. 
Its  value  is  derived  from  his  own  work  (Park  (1990))  concerning  the  partition  function  of 
the  involved  species.  In  addition,  all  radiative  processes  have  been  removed  in  our  CR 
model  in  order  to  test  the  kinetic  scheme  only. 

The  situations  adopted  for  illustrating  the  difference  between  these  schemes  and  our  CR 
model  are  related  to  two  pressure  levels  of  2  kPa  and  10  kPa  framing  the  typical  range 
over  which  a  re-entry  occurs.  I11  both  cases,  the  plasma  is  initially  at  equilibrium  with 
Te  =  Ta  =  7000  I\  and  is  recombining  at  Te  =  Ta  =  3000  A'.  We  can  see  on  Fig.  20  the 
time  evolution  of  the  neutral  species  and  that  of  the  charged  species  except  electrons  on 
Fig.21  calculated  with  our  CR  model  and  with  the  three  other  kinetic  schemes.  Figures 
22  and  23  present  the  calculated  evolutions  in  the  same  conditions  but  for  p  —  10  kPa. 
Figure24  illustrates  finally  the  time  evolution  of  the  electron  density  in  the  same  thermal 
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conditions  as  previously  considered  with  p  —  2  kPa  (a)  and  with  p  —  10  kPa  (/ 3 ). 


Figure  20:  Time  evolution  of  JV2,  02,  NO, 
N  and  O  in  a  typical  re-entry  situation:  re¬ 
combination  at  p  =  2  kPa  from  Te  =  T 4  = 
7000  A'  to  Te  =  Ta  =  3000  A”  (thick  solid 
line:  our  model,  thin  solid  line:  Dunn  & 
Kang,  thin  dotted  line:  Gupta,  thick  dot- 
dashed  line:  Park). 


Figure  21:  Time  evolution  of  N^,  0£, 
N 0+ ,  N+  and  0+  in  a  typical  re-entry  sit¬ 
uation:  recombination  at  p  =  2  kPa  from 
Te  =  Ta  =  7000  K  to  Te  =  TA  =  3000  K 
(thick  solid  line:  our  model,  thin  solid  line: 
Dunn  &  Kang,  thin  dotted  line:  Gupta, 
thick  dot-dashed  line:  Park). 


Since  the  plasma  considered  in  this  section  is  optically  thick  in  order  to  make  a  compar¬ 
ison  between  chemical  processes  involved  by  the  kinetic  schemes  only,  the  final  state  is 
an  equilibrium  state.  The  agreement  is  rather  good  for  this  final  equilibrium  state  con¬ 
cerning  N2,  O2,  NO,  N  and  O  except  for  the  simplified  kinetic  scheme  of  Gupta  which 
leads  to  a  final  overpopulation  of  N  and  to  a  final  underpopulation  of  N2  and  NO.  A 
calculation  under  a  higher  pressure  condition  as  the  atmospheric  one  show  the  same  trend 
with  in  addition  an  overpopulation  of  02-  This  behavior  is  the  direct  result  of  the  influ¬ 
ence  of  the  equilibrium  constants.  In  particular,  those  of  the  Zeldovich’s  reaction  (32a)  is 
lower  leading  to  a  lower  value  for  the  rate  coefficient  in  the  case  of  the  kinetic  scheme  of 
Gupta:  the  population  density  of  02  and  N  are  then  higher  while  those  of  iV2  and  NO 
are  lower.  We  can  observe  for  the  ions  on  Figures  21  and  23  that  their  population  density 
is  of  the  same  order  of  magnitude  except  for  O 2  and  NO+ .  At  atmospheric  pressure, 
02  is  overpopulated  in  the  final  equilibrium  state  in  the  case  of  the  kinetic  scheme  of 
Gupta.  This  is  also  the  direct  result  of  the  different  equilibrium  constants  involved  by 
each  kinetic  scheme.  Moreover,  we  can  observe  that  NO+  is  the  more  significant  ion  in 
the  final  equilibrium  state:  the  temperature  being  relatively  low,  the  number  density  of 
NO+  is  the  highest  one  as  a  result  of  the  high  concentration  of  NO  and  its  relatively  low 
ionization  characteristic  temperature  (~  107500  A'). 

The  final  electron  density  is  in  strong  disagreement  from  a  kinetic  scheme  to  another. 
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Figure  22:  Time  evolution  of  JV2,  O2,  NO, 
N  and  O  in  a  typical  re-entry  situation:  re¬ 
combination  at  p  =  10  kPa  from  Te  =  TA  = 
7000  K  to  Te  =  Ta  =  3000  K  (thick  solid 
line:  our  model,  thin  solid  line:  Dunn  & 
Kang,  thin  dotted  line:  Gupta,  thick  dot- 
dashed  line:  Park). 


Figure  23:  Time  evolution  of  1V^,  O^, 
N 0+ ,  N+  and  0+  in  a  typical  re-entry  sit¬ 
uation:  recombination  at  p  =  10  kPa  from 
Te  =  TA  =  7000  K  to  Te  =  TA  =  3000  K 
(thick  solid  line:  our  model,  thin  solid  line: 
Dunn  &  Kang,  thin  dotted  line:  Gupta, 
thick  dot-dashed  line:  Park). 


Figure  24:  Time  evolution  of  the  electron  density  for  a  recombining  plasma  from  Te  =  TA  = 
7000  K  to  Te  =  Ta  =  3000  K  with  p  =  2  kPa  (a)  and  with  p  =  10  kPa  (/ 3 )  (thick  solid  line:  our 
model,  thin  solid  line:  Dunn  &:  Kang,  thin  dotted  line:  Gupta,  thick  dot-dashed  line:  Park). 
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The  reason  can  be  ascribed  once  again  to  the  equilibrium  constants,  but  also  to  the  elec¬ 
troneutrality:  since  NO+  has  the  highest  number  density  among  the  ions  and  since  our 
kinetic  scheme  provides  the  highest  population  density  for  this  molecular  ion,  the  value 
of  ne  in  the  final  equilibrium  state  is  then  the  highest  one. 

Although  the  final  state  is  related  with  too  long  time  scales  with  respect  to  the  conditions 
typically  involved  in  the  re-entry  conditions,  nevertheless  its  study  is  very  important:  this 
informs  about  the  quality  of  the  equilibrium  constants  which  are  crucial  because  they 
provide  the  calculation  of  the  inverse  rate  coefficients  and  consequently  they  influence 
the  time  evolution  of  the  species  number  densities.  The  comparison  between  the  good 
correlation  of  the  final  state  for  the  neutral  species  and  the  bad  one  for  that  related  to  the 
charged  species  lead  us  to  conclude  that  the  equilibrium  constants  for  the  neutral  species 
are  better  than  for  the  charged  ones. 

Except  the  influence  of  the  rate  coefficients,  the  previous  conclusion  explain  in  part  the 
fact  that  the  different  population  densities  evolve  strongly  in  a  different  way  during  the 
transient  step.  We  have  already  mentioned  the  influence  of  the  rate  coefficient  for  the 
Zeldovich’s  reaction  (32a)  which  is  very  different  from  a  kinetic  scheme  to  another.  The 
equilibrium  constant  being  also  different  from  a  kinetic  scheme  to  another,  it  is  clear  that 
the  time  evolution  has  no  reason  to  be  the  same  one.  Nevertheless,  the  discrepancy  is 
sometimes  in  a  ratio  higher  than  50.  For  relatively  short  time  scales  as  those  related  to 
typical  re-entry  (~  1CR4s),  the  concentrations  of  O2  and  NO  are  very  different. 

In  addition,  this  evolution  depends  strongly  on  the  adopted  initial  conditions:  equilibrium 
at  T  =  7000  K  induces  here  higher  concentrations  for  N+  and  0+  than  for  NO+ :  as  a 
result,  the  latter  ion  plays  a  relatively  minor  role.  For  a  lower  temperature  as  considered 
in  the  next  section,  one  have  [ NO+ }  =  4.9  x  1018m-3  at  Tc  =  Ta  =  5400  K  while  [ N+ ]  = 
1.9  x  1017m-3  and  [0+]  =  2.4  x  1017m-3  when  p  =  2  kPa  and  [NO+]  =  1.9  x  1019m-3, 
[N+]  =  1.4  x  1017m-3  and  [0+]  =  3.5  x  1017m-3  when  p  =  10  kPa.  In  this  case,  NO+ 
is  the  more  important  ion  and  the  recombination  of  electrons  is  driven  essentially  by  its 
DR. 

Figure  24  illustrates  the  evolution  of  ne :  the  different  kinetic  schemes  are  globally  in 
strong  disagreement,  even  if  the  schemes  of  Dunn  &  Kang  and  Gupta  are  relatively  well 
correlated  for  relatively  short  time  scales.  They  underestimate  clearly  the  electron  density 
in  relation  to  that  calculated  by  the  other  schemes  for  t  <  ICR4  s  which  is  of  interest  in 
the  case  of  a  typical  re-entry.  In  the  case  electrons  play  a  role  the  temperature  being 
sufficiently  high,  the  kinetic  scheme  chosen  to  solve  the  boundary  layer’s  equations  will 
lead  to  very  different  results.  We  have  made  a  special  effort  to  describe  correctly  the 
recombination  and  as  a  result  the  ionization  via  the  equilibrium  constants:  we  can  affirm 
that  the  recombination  is  too  strong  with  the  three  common  simplified  kinetic  schemes 
here  studied. 

Finally,  Figure  19  shows  that  the  relaxation  of  the  vibrational  distribution  for  A^O+(X4S+) 
is  almost  complete  when  the  transient  regime  begins  for  the  other  species:  the  assump¬ 
tion  of  a  Maxwell  equilibrium  between  the  vibration  of  molecules  with  the  other  species 
is  validated. 
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4.4  Conclusions 

Three  main  conclusions  can  be  derived  from  the  study  of  the  collisional  radiative  model 
for  air.  First,  ionization  situations  have  shown  that  and  0~ ,  due  to  their  unsteady 
behavior,  may  have  an  influence  on  the  characteristic  time  scale  necessary  to  reach  the 
equilibrium  state.  The  common  approach  is  to  ignore  these  ions  in  shock  layer  computa¬ 
tions.  ft  could  be  interesting  to  account  for  them  in  a  shock  calculation. 

The  second  important  conclusion  is  that  a  CR  model  has  to  be  time-dependent:  it  provides 
chemical  characteristic  time  scales  which  can  be  compared  with  the  other  characteristic 
time  scales  of  the  flow.  This  temporal  approach  is  also  particularly  powerful  to  check  the 
quality  of  the  equilibrium  constants  (by  the  equilibrium  state  finally  reached  and  also  by 
the  transient  regime). 

Third,  the  comparison  between  the  three  common  simplified  kinetic  schemes  used  for 
re-entry  studies  and  the  detailed  CR  model  has  shown  that  the  results  are  very  different 
from  a  scheme  to  another.  This  is  the  reason  why  it  is  clear  that  the  better  solution  to 
validate  a  kinetic  scheme  is  in  fact  to  study  the  transient  regime.  We  have  shown  that 
the  final  state  allows  to  test  the  quality  of  the  equilibrium  constants  only.  The  next  step 
of  this  work  is  to  derive  a  new  simplified  kinetic  scheme  to  be  implemented  in  a  Navier- 
Stokes  code  devoted  to  re-entry  studies.  A  time-dependent  approach  will  allow  to  verify 
the  relevance  of  this  new  reduced  kinetic  scheme  over  a  wide  temporal  range. 


5  Stagnation  point  boundary  layer 

A  method  has  been  proposed  by  researchers  (Kolesnikov  (1999))  at  IPM  (Moscow)  to 
determine  the  catalytic  properties  of  thermal  protection  systems  (TPS).  This  method  is 
indirect  and  is  based  on  comparisons  of  the  wall  heat  flux  measured  on  a  TPS  sample 
set  in  a  high  enthalpy  flow  generated  by  a  ICP  torch  with  the  flux  calculated  using  a 
stagnation  point  boundary  layer  code.  This  method  is  popular  and  is  used  in  different 
places  in  Europe,  for  example,  at  CORIA  (France),  VKI  (Belgium)  and  in  EADS-ST 
(France).  Studies  have  been  carried  out  for  different  materials.  Most  results  available 
in  the  literature  are  for  pure  silica  or  for  material  containing  mostly  silica  (e.g.  RCG 
contains  92%  silica).  In  this  section,  we  present  the  stagnation  point  boundary  layer  code 
for  air  with  different  models  for  catalycity  effects.  Then  for  two  typical  flow  conditions 
encountered  in  catalycity  studies,  the  influence  of  the  catalycity  model  and  of  the  kinetic 
scheme  is  disccused. 


5.1  Stagnation  point  boundary  layer  code 

Classical  boundary  layer  equations  are  written  assuming  that  the  flow  is  axisymmetric 
and  laminar,  dissociated  but  non  ionized,  in  chemical  nonequilibrium  and  in  thermal 
equilibrium.  Using  the  classical  Lees-Dorodnitsyn  (x,  y )  — >  (£,  77)  transformation,  the 
boundary  layer  equations  at  the  stagnation  point  reduce  to  a  set  of  ordinary  differential 
equations  (Barbante  (2001)) 


•  continuity 


dV 

dr\ 


+  f  =0 


(50) 
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momentum 


species  i  continuity 


df_  =  1  (ps  _  ,2)  d(jdf_ 
dp  2  \  p  J  dr]  \ 0  dr] 


v^  +  aS  =  w, 


dr]  dr] 


(51) 


(52) 


energy 


ydg_  _  d_  /  l^dg_ 


with 


and 


9  = 


dyi  hi  \ 
Pr  dp  hs  J 

d 

dp 

(53) 

u 

(54) 

Us 

h 

hs 

(55) 

where  the  subscript  5  indicates  a  quantity  evaluated  at  the  boundary  layer  outer  edge,  yi 
is  species  i  mass  fraction,  u  is  the  velocity  component  parallel  to  the  wall  and  h  is  the 
enthalpy  of  the  mixture.  Iq  is  defined  as 

Ph 


In  — 


PsPs 


The  local  Prandtl  number  is 

The  dimensionless  diffusion  flux  is 


Pr  = 


J?  = 


j-LCp 

— 

J? 


2 


(56) 

(57) 

(58) 


and  the  dimensionless  mass  production  term  is 

m=  w‘ 


2 


(59) 


where  Wi  is  the  mass  source  term  of  species  i  in  the  (x,y)  framework. 

In  order  to  carry  out  a  stagnation  point  boundary  layer  computation,  it  is  necessary  to 
know  the  wall  temperature  Tw  and  three  parameters  at  the  boundary  layer  edge:  the 
velocity  gradient  (3  =  T$  and  Ps. 

It  is  important  to  mention  that  to  do  comparisons  with  catalycity  studies  in  subsonic  flows 
generated  by  an  ICP  torch,  the  finite  thickness  of  the  boundary  layer  at  low  Reynolds 
numbers  has  to  be  taken  into  account  (Barbante  (2001)).  Then,  the  set  of  equations 
for  the  stagnation  point  boundary  layer  has  to  be  slightly  modified.  These  corrections 
require  to  know  precisely  the  boundary  layer  thickness  and  different  velocity  gradients. 
In  this  work,  we  propose  to  discuss  the  influence  of  different  models  on  the  results  of  the 
stagnation  point  boundary  layer  simulation,  and  therefore  these  corrections  have  been 
neglected. 
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5.1.1  Catalycity  and  boundary  conditions 

Catalycity  effects  are  generally  taken  into  account  using  two  types  of  dimensionless  co¬ 
efficients:  a  recombination  probability  7  (0  <  7  <  1)  and  an  energy  accommodation 
coefficient  at  the  wall  (3.  The  7*  coefficient  for  a  species  i  is  defined  as 


(60) 


where  Fi  is  the  flux  of  species  i  impinging  on  the  wall,  and  Fi^recomb  is  the  recombining 
flux  of  species  i  on  the  wall.  Assuming  a  Boltzmann  distribution  of  velocities  close  to  the 

wall 


(61) 


where  m,  is  the  mass  of  species  i.  Then,  the  production  rate  of  species  i  at  the  wall  Wi}Cat 
can  be  written  as 


(62) 


The  (3  coefficient  is  defined  as 


(63) 


where  q  is  the  effective  energy  flux  given  to  the  wall,  D  is  the  dissociation  energy  of  the 
molecule  formed.  Then  the  heat  flux  transferred  to  the  wall  depends  on  the  product 


7*  =  7 f3.  It  is  interesting  to  note  that  7'  and  7  are  equal  when  the  time  of  stay  of  the 


formed  molecules  on  the  surface  is  long  enough.  I11  this  case,  the  whole  recombination 
energy  is  given  to  the  wall  (/3  =  1).  In  fact,  there  has  been  only  a  few  studies  on  the  (3 
coefficient  determination  and  its  modeling  (Daiss  et  al.  (1997);  Scott  (1990);  Suslov  and 
Tirskyi  (1994);  Barbato  et  al.  (1998)),  therefore  generally  this  coefficient  is  assumed  to 
be  equal  to  1  and  only  the  7  coefficient  is  used.  In  this  work,  we  have  also  used  this 
approximation.  Then,  boundary  conditions  are  written  as 
At  the  wall  (77  =  0) 


(64) 

(65) 

(66) 


u  —  v  —  0 


where  Jfw  is  the  diffusion  flux  of  species  i  in  the  direction  normal  to  the  wall.  At  the 
boundary  layer  edge  (77  — >  00) 


(67) 

(68) 
(69) 


U  =  Us 

h  =  hs  =  h(T§,  Vi  s) 

Vi  =  Vi, s 


The  boundary  layer  flow  is  assumed  to  be  in  chemical  equilibrium  at  the  boundary  layer 
edge,  then  yits  =  yi(T& ,  P$).  Finally,  the  wall  heat  flux  is  given  by 


W 


(70) 
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5.1.2  Thermochemical  model 

In  the  stagnation  point  boundary  layer  code,  air  is  considered  as  a  5  species  mixture  (N2, 
O2,  NO,  N  and  O ).  To  determine  the  viscosity  and  the  thermal  conductivity  of  the 
mixture,  the  mixing  laws  given  by  Gupta  et  al.  (1990)  have  been  used  with  the  collision 
integrals  calculated  recently  by  Capitelli  et  al.  (2000b).  For  boundary  layer  flows,  it  is 
particularly  important  to  calculate  accurately  diffusion  fluxes.  The  classical  Fick’s  law 
gives  a  simple  expression  of  diffusion  fluxes.  For  example,  the  y  component  of  the  diffusion 
flux  of  species  i  in  the  mixture  is 

J?  =  (71) 

where  Dim  is  the  diffusion  coefficient  of  species  i  in  the  mixture.  However,  this  simple 
expression  does  not  satisfy  the  mass  conservation:  J?  =  0  (unless  all  the  Dirn 

coefficients  are  equal)  and  therefore  should  not  be  used  when  diffusion  processes  are 
important.  The  rigorous  expression  of  diffusion  fluxes  is  given  by  Stefan  -  Maxwell  and 
Sutton  and  Gnoffo  (1998)  have  proposed  an  efficient  iterative  technique  to  determine 
them.  In  this  case,  the  calculated  fluxes  satisfy  the  constraint  of  mass  conservation. 

In  this  work,  the  kinetic  schemes  proposed  by  Dunn  and  Kang  (1973),  Park  (1993)  and 
Gupta  et  al.  (1990)  adapted  to  a  five  species  air  mixture  have  been  used. 


5.1.3  Numerical  method 


The  boundary  layer  equations  (except  the  continuity  equation)  can  be  written  in  the 
generic  form 


with  boundary  conditions 


82f  r.9f  r  j 

+  c/  =  d 

orjz  or 1 

(72) 

n  df  f 

77  =  0  :  ai—  +  bif  =  ci 

(73) 

V  =  Vs  :  /  =  fs 

(74) 

and  are  solved  with  a  4th  order  finite  difference  method  (Barbante  (2001)).  500  points 
are  taken  in  the  rj  direction. 


5.2  Catalycity  models 

5.2.1  Simplified  models 

The  simplest  model  is  to  consider  that  only  two  recombination  reactions  occur  at  the 
wall  (i.e.  N  +  N  — >  N2  and  O  +  O  — >  O 2)  and  that  these  two  reactions  have  the 
same  7  coefficient.  This  model  is  widely  used  to  determine  an  average  7  coefficient  for  a 
material.  From  different  measurements  in  pure  nitrogen  and  oxygen  flows,  Scott  (1990) 
and  Kolodziej  and  Stewart  (1987)  have  determined  recombination  rates  of  N  in  N2  and 
O  in  02  on  different  surfaces  as  a  function  of  the  wall  temperature  Tw.  Scott’s  fits  for 
recombination  coefficients  are 

7 NN  =  0.0714e-2219/T”,  950  <  Tw  <  1670/1  (75) 
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700  =  16.e-10271/T”,  1400  <TW  <  1650 K  (76) 

Kolodziej  and  Stewart  (1987)  give 

7 NN  =  0.061e-2480/T“,  1410  <TW  <  1640/1  (77) 

700  =  40.e-11440/T“,  1435  <TW  <  1580 K  (78) 

These  rates  are  usually  implemented  in  air  flow  codes  assuming  that  no  NO  is  formed  on 
the  surface.  They  are  very  simple  to  use,  but  are  valid  in  a  small  temperature  range. 

5.2.2  Phenomenological  models 

Different  authors  (Kovalev  and  Suslov  (1996),  Jumper  and  Seward  (1994),  Nasuti  et  al. 
(1996),  Deutschmann  et  al.  (1995)  and  Daiss  et  al.  (1997))  have  developed  phenomeno¬ 
logical  models  in  which  the  gas-surface  interaction  is  represented  by  different  elementary 
reactions 

•  atom  adsorption/desorption  (Fig.  25) 

•  recombination  of  an  atom  of  the  gas  with  an  atom  adsorbed  on  the  wall:  Eley-Rideal 
(E-R)  type  process  (Fig.  26) 

•  recombination  of  two  adsorbed  atoms  at  the  wall:  Langmuir-Hinshclwood  (L-H) 
type  process  (Fig.  27) 

For  the  temperatures  and  pressures  studied,  molecule  adsorption  at  the  wall  is  neglected. 
A  rate  coefficient  is  associated  to  each  elementary  process.  In  this  case,  the  global  7 
coefficient  depends  on  temperature  but  also  on  the  gas  composition  close  to  the  wall. 


A 


Figure  25:  atom  adsorption/desorption 


Figure  26:  Eley-Rideal  (E-R)  type  recom¬ 
bination  process 


Figure  27:  Langmuir-Hinshelwood  (L-H) 
type  recombination  process 
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In  the  model  derived  by  Deutschmann  et  ah,  simplified  expressions  for  elementary  rate 
coefficients  are  used.  In  works  of  Jumper  et  al.,  Daiss  et  ah  and  Nasuti  et  al.,  a  better 
description  of  each  process  is  given.  In  the  model  of  Jumper  et  ah,  recombinations  of  O 
in  02  and  N  in  Ah  are  assumed  to  be  Eley-Rideal  type  processes.  In  the  models  of  Daiss 
et  al.  and  Nasuti  et  ah,  recombination  reactions  are  assumed  to  be  of  both  types.  It  is 
interesting  to  note  that  the  model  of  Nasuti  et  ah  takes  into  account  the  NO  formation 
at  the  wall.  In  the  model  of  Nasuti  et  al,  the  following  processes  are  taken  into  account 


•  N  and  O  atom  adsorptions  and  desorptions 

O  +  s^O  —  s,  N  +  s^N  —  s  (79) 

•  Eley-Rideal  type  processes 

O  —  s  +  0^02  +  s,  N  —  s  +  N^N2  +  s  (80) 

0-s  +  N<->  NO +  s,N-s  +  0<->  NO +  s  (81) 

•  Langmuir-Hinshelwood  type  processes 

O  —  s  +  O  —  s  ^  02  +  2s,  N  —  s  +  N  —  s^>N2  + 2s  (82) 

0-s  +  N-s^N0  +  2s  (83) 


where  s  is  a  free  surface  site.  With  this  approach,  it  is  necessary  to  define  four  recombi¬ 
nation  probabilities  at  the  wall  :  7 OG,  ^nnNnoiIon- 

Expressions  of  7  coefficients 

Due  to  a  competition  between  O  and  N  recombination  processes  at  the  wall,  the  four 
recombination  probabilities  are  not  independent  of  each  other.  Therefore  7 nn  and  7 00 
coefficients  obtained  in  pure  nitrogen  and  oxygen  flows  certainly  overestimate  these  coef¬ 
ficients  in  an  air  flow  where  both  atoms  are  present.  Furthermore,  the  number  of  surface 
sites  for  N  and  O  adsorption  may  be  limited  due  to  the  chemisorption  on  the  surface  of 
other  species  such  as  the  H20  molecule  which  forms  Si-OH  bounds  on  a  silica  surface.  To 
determine,  the  surface  coverage  rate  by  OH ,  Nasuti  et  al.  proposes 

90h  =  -4(1  -  exp{-EOH/kBT ))  (84) 

The  surface  coverage  rates  0O  et  9N  are  related  to  the  free  surface  sites  9f  by 

9f  =  1  —  9o  —  9n  —  9oh  (85) 

with  0  <  9f  <  1.  At  the  stationary  state,  the  number  of  adsorbed  atoms  on  the  surface  is 
constant  and  therefore,  the  flux  of  adsorbed  atoms  i  on  the  surface  ( Fai )  is  equal  to  the 
flux  of  atoms  recombining  following  a  L-H  ( Fihii )  or  a  E-R  ( Ferii )  process  and  to  the  flux 
of  atoms  thermally  desorbed  (Ftdii).  Then,  for  O  and  N  we  have 

FaO  =  Fer]\fo  +  Feroo  +  FihNO  +  FihOO  +  Ftdo  (86) 
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FaN  —  FerON  +  FerNN  +  Fit iNO  +  FihNN  +  FtdN  (87) 

As  we  will  see  in  the  following,  these  different  fluxes  are  functions  of  the  two  variables 
0o  and  dN.  The  system  of  equations  (86)  and  (87)  is  therefore  a  quadratic  system  of 
equations  for  Oq  and  0 n  variables  which  can  be  easily  solved  by  an  iterative  method 
(e.g.  Newton- Raphson  method).  When  Oo  et  On  are  known,  the  different  recombination 
probabilities  can  be  expressed  as 


7  oo  = 


Inn  — 


( Fihoo  +  2  Feroo) 
F~o 

( F^nn  +  2FerNN) 


F 


N 


Ion  — 


(FihNO  +  FerNO  +  FerQN) 


Ft 


o 


iNO  ~ 


( F^no  +  FerNO  +  Fer0N ) 


N 


These  four  coefficients  depend  on  pQl  Pn ,  On,  0o  and  T. 


(88) 

(89) 

(90) 

(91) 


Wall  boundary  conditions  for  the  different  species 

The  production  rate  WitCat  at  the  wall  of  species  i  =  N2,  02,  NO,  N  and  O  can  be  written 
as 


Wo  =  m0(pioo  +  1on)F0 

(92) 

Wn  =  ttin{inn  +  1no)Fn 

(93) 

wo2  =  -mo^ooFo 

(94) 

=  —(mo'yoNFo  +  itin^noFn) 

(95) 

wn2  —  ~m nInnFn 

(96) 

where  Fi  is  the  flux  of  species  i  impinging  on  the  wall  and  is  given  by  Eq.  (61).  The  flux 
of  atoms  i  adsorbed  at  the  wall  Fai  is  proportional  to  Ft  and  to  the  number  of  free  surface 
sites  Of 

Fai  =  SiOfFi  (97) 


where  s*  is  a  sticking  coefficient  (0  <  s*  <  1).  The  flux  Ferij  of  atoms  i  in  the  gas, 
recombining  with  adsorbed  atoms  j  to  form  an  ij  molecule  following  an  E-R  process  is 
proportional  to  E)  and  0j 


where  0  <  7*-  <  1  with: 

'  LJ 


F  =  ^*0F- 

1  enj  hjwJ±  7 

r!ij  Perij&^Cpi,  Qerij/^BP) 


(98) 

(99) 


In  a  L-H  recombination  process,  both  atoms  are  adsorbed  at  the  wall.  The  flux  of  ad¬ 
sorbed  atoms  which  recombine  following  a  L-H  type  process  to  form  02,  N2  and  NO  are 
respectively 

Zmoo  =  2  von02oexp[-{QihoolkBT)\  (100) 

ZihNN  =  2uNn0^exp[— {Qihnn /^bT)]  (101) 
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ZihNo  —  (ov  +  iyo)n^N^oexP[^  (QihNo/kBT)\  (102) 

where  n  is  the  number  of  surface  sites  per  unit  surface  and  z/,;  is  a  characteristic  frequency 
given  by 

V%  =  (cA/ A)  a/7t  kBT  /  2m,i  (103) 

where  A  is  the  mean  distance  between  surface  sites  and  cA  is  a  constant.  Activation 
energies  are  given  by 

Qihoo  —  M AX(2Qa0  —  D0ch  Em0)  (104) 

where  Qai  is  the  surface-atom  bond  energy  for  atom  i  and  Erm  is  the  energy  necessary  for 
species  i  to  move  on  the  surface  and  Dij  is  the  bond  energy  between  i  and  j  atoms. 


QihNN  —  MAX(2Q , 

in  ~  Dnn,  Em N) 

(105) 

QlhNO  ~  QaO  +  QaN  ~  D^o 

The  flux  of  thermal  desorption  Ftdi  is 

(106) 

Ftdi  =  nOi(kBT/h)exp[-Qtdi/kBT )] 

(107) 

where  Qt(a  is  an 

activation  energy. 

For  N 

SN 

PerNN 

QerNN 

QtdN 

0.015 

0.1 

20  kJ/mol 

QaN 

Wall 

n 

A 

4.51018m"2 

510 10m 

For  O 

For  N 

so 

P °rOO 
QerOO 
QtdO 

0.05 

0.1 

29k  J/mol 

Qao/2. 

EmN 

QaN 

Dnn 

82.2kJ  /mol 

530.8  kJ /mol 

94:1.3k  J/mol 

For  0 

CA 

3.5 

EmO 

QaO 

Dqo 

159.1  kJ/mol 

499.8 k  J/mol 

498 kJ/mol 

For  OH 

A 

Eqh 

1. 

2.5  kJ/mol 

Table  10:  Data  for  the  model  derived  by 
Nasuti  et  al.  (1996). 

Table  11:  Parameters  derived  to  have 
the  best  agreement  with  experimental 
results  obtained  in  pure  nitrogen  and 

oxygen  flows  (Nasuti  et  al.  (1996)). 


Determination  of  the  parameters  of  the  model 
Among  the  different  parameters  of  the  model,  there  are: 

•  parameters  which  depend  on  the  model  chosen  for  the  surface  :  Qai,  A,  n  and  Emi 
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•  fundamental  data  available  in  the  literature  as  bond  energies  Dl3 

•  parameters  which  can  be  derived  from  semi-empirical  models  :  Qerij,  Qai  and  Qtdi, 

•  parameters  for  which  no  theoretical  model  is  available:  Perij-,  ,7  and  <74. 

The  three  first  types  of  parameters  are  given  in  table  10.  Figure  28  shows  numerous 
measurements  reviewed  by  Jumper  and  Seward  (1994)  and  Wiley  (1993)  of  7 nn  et  7 00 
for  pure  silica  or  RCG  in  pure  nitrogen  and  oxygen  flows.  These  results  have  been  used 
to  determine  the  parameters  of  the  model  given  in  table  11.  For  the  thermal  desorption 
threshold,  we  have  assumed  that  QtdN  =  QaN  and  Qtd,o  is  taken  equal  to  Qao/ 2  to  have 
a  maximum  of  7 00  at  T  ~  1600/1.  On  Fig.  28,  the  results  obtained  by  the  model 
proposed  by  Jumper  are  also  shown.  In  this  model,  recombination  processes  are  of  the 
E-R  type  and  Perij  and  s*  are  assumed  to  depend  on  the  temperature.  Experimental 
results  seem  to  indicate  that  7 00  has  a  maximum  at  T  ~  1600/1.  This  evolution  is 
well  reproduced  by  both  models  and  is  due  to  thermal  desorption  which  becomes  very 
efficient  as  temperature  increases.  With  the  model  of  Nasuti,  7 00  is  independent  of 
pressure  for  T  <  lOOO/i  and  depends  on  the  pressure  for  higher  temperatures.  For 
O  in  O2  recombination,  we  have  checked  that  in  the  studied  wall  temperature  range 
300 K  —  2000/1 ,  the  L-H  type  process  is  negligible,  and  then  the  recombination  process  is 
an  E-R  type  recombination  process.  For  'Jnn,  experimental  data  are  scarce  and  therefore 
the  existence  of  a  maximum  at  T  ~  1600/1  is  not  obvious.  This  maximum  exists  in  the 
model  of  Jumper  but  not  in  the  model  of  Nasuti.  With  the  model  of  Nasuti,  depends 
on  pressure  only  for  600 K  <  Tw  <  1000/1 .  For  T  <  600 K,  the  recombination  process  is 
an  E-R  type  process  and  it  is  a  L-H  process  for  T  >  1000/1 .  For  air,  no  measurements 
of  the  different  7  coefficients  are  available.  For  the  recombination  of  O  in  02  and  of  N  in 
N'2,  the  parameters  derived  hereabove  are  used.  The  other  parameters  are  given  in  table 
12  following  the  assumptions  of  Nassuti  et  al. 


Figure  28:  Recombination  probabilities  of  N  in  1V2  and  of  O  in  02  on  pure  silica  and  RCG. 
Symbols  correspond  to  different  experimental  results  reviewed  in  Jumper  and  Seward  (1994) 
and  Wiley  (1993). 
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For  ON 

QerON 

QerNN 

PerON 

PerNN 

For  NO 

QerNO 

QerOO 

Per  NO 

PerOO 

Dno 

625.1  kJ/mol 

Table  12:  Parameters  for  ON  and  NO  (Nasuti  et  al.  (1996)). 


5.3  Results 

In  this  section,  we  present  results  for  two  typical  conditions  of  catalycity  studies  carried 
out  in  high  enthalpy  flows  generated  by  a  ICP  torch.  For  both  test-cases,  we  propose  to 
discuss  the  influence  of  the  model  used  for  catalycity  and  for  the  gas  kinetic  scheme  on 
the  calculated  results. 

5.3.1  High  pressure  test  case  :  P  =  lOOOOPa 

We  consider  typical  boundary  edge  conditions  at  P$  =  lOOOOPa:  T$  =  4500 K  and 
f3  =  11600s-1  (Barbante  (2001)).  The  simplest  model  for  catalycity  presented  in  Sec¬ 
tion  5.2.1  is  to  consider  that  only  two  recombination  reactions  occur  at  the  wall  (i.e. 
N  +  N  — >  N2  and  O  +  O  ^  O2)  and  that  these  two  reactions  have  the  same  7  coefficient. 
Figure  29  presents  the  wall  heat  flux  as  a  function  of  the  wall  temperature  for  different 
values  of  7.  In  the  classical  method  derived  by  IPM  to  determine  the  catalycity  of  a  TPS, 
measured  heat  fluxes  are  set  in  this  figure  to  determine  the  7  coefficient  of  the  studied 
wall. 

As  expected  the  calculated  wall  heat  flux  increases  as  7  increases  and  as  the  wall  temper¬ 
ature  decreases.  For  a  given  wall  temperature,  the  heat  flux  in  the  fully  catalytic  case  is 
about  twice  the  flux  in  the  noncatalytic  case.  The  heat  flux  for  a  fully  frozen  flow  (  i.e. 
no  reactions  in  the  gas  and  at  the  wall)  is  also  shown.  This  flux  is  the  lower  heat  flux 
limit.  In  the  noncatalytic  case,  the  difference  between  the  flux  calculated  in  the  frozen 
case  and  the  curve  denoted  7  =  0,  shows  the  influence  of  the  gas  chemistry  on  the  heat 
flux.  Figure  29  shows  that  this  influence  decreases  as  the  wall  temperature  increases, 
as  the  recombination  of  atoms  in  the  gas  is  more  efficient  at  low  temperature.  Figure 
30  presents  the  wall  heat  flux  calculated  using  the  phenomenological  model  derived  by 
Nasuti  et  al  and  presented  in  Section  5.2.2  .  At  low  temperature,  this  model  shows  that 
the  wall  is  noncatalytic  and  that  its  catalycity  increases  up  to  Tw  ~  1400/1 .  For  higher 
temperatures,  thermal  desorption  becomes  very  efficient  and  the  calculated  wall  heat  flux 
decreases. 
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Figure  29:  Wall  heat  flux  as  a  function 
of  the  wall  temperature  for  different  values 
of  the  average  7  coefficient.  The  kinetic 
scheme  of  Dunn  and  Kang  is  used. 


Figure  30:  Same  as  Fig.  29  with  the  heat 
flux  calculated  using  the  phenomenological 
model  derived  by  Nasuti  et  al. 


Figure  31  shows  that  the  7  coefficients  of  the  different  species  calculated  with  the  phe¬ 
nomenological  model  have  very  different  values  from  each  other  and  different  evolutions 
with  the  wall  temperature.  At  low  temperature  Tw  <  1000/1 ,  7 00  and  7jvo  are  close  to 
each  other  and  much  higher  than  the  two  other  coefficients.  The  coefficient  increases 
as  the  wall  temperature  increases  and  becomes  the  highest  recombination  coefficient  for 
Tw  >  1800/1 .  We  note  that  the  70 n  remains  very  small  in  the  whole  temperature  range. 


Figure  31:  7  coefficients  of  the  phe¬ 

nomenological  model  derived  by  Nasuti  et 
al.  for  the  conditions  of  Fig.  30. 


Figure  32:  Different  values  of  0  coefficients 
for  the  same  condition  as  Fig.  31. 


Figure  32  shows  the  surface  coverage  for  O,  N  atoms  and  the  OH  molecule  for  the 
conditions  of  Fig.  31.  The  surface  coverage  by  the  OH  molecule  decreases  as  the  wall 
temperature  increases.  The  surface  coverage  by  O  atoms  is  on  the  average  of  0.5  for 
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300/1  <TW<  1500/1  and  then  decreases  rapidly  for  higher  temperature  due  to  thermal 
desorption.  The  N  atom  surface  coverage  remains  very  small  in  the  whole  temperature 
range. 

Figure  33  shows  the  species  concentration  for  Tw  =  1500/1  and  the  phenomenological 
model  proposed  by  Nasuti  et  al.  Strong  variations  are  observed  close  to  the  wall  for  N, 
O2  and  NO.  Figure  33  compares  the  wall  heat  flux  calculated  with  the  phenomenological 
model  proposed  by  Nasuti  et  al.  using  different  kinetic  schemes  (Dunn  and  Kang  ,  Park 
and  Gupta)  and  for  a  frozen  flow.  In  the  scheme  of  Park,  backward  rate  coefficients  are 
derived  from  forward  rates  using  equilibrium  constants  fitted  as  a  function  of  temperature 
for  T  >  500/1 .  To  avoid  numerical  problems  at  low  temperatures  with  these  fits,  we 
present  results  for  the  kinetic  scheme  of  Park  for  T  >  1000/1 .  The  heat  fluxes  calculated 
using  the  kinetic  scheme  of  Dunn  and  Kang  are  very  close  to  the  results  obtained  for 
the  scheme  of  Gupta  for  T  >  700/1 .  The  most  significant  discrepancies  are  observed  for 
lower  temperatures.  The  heat  fluxes  calculated  using  Park’s  kinetic  scheme  are  slighlty 
higher  than  those  obtained  with  Dunn  and  Kang  scheme  but  the  discrepancy  remains 
small.  The  heat  flux  calculated  assuming  a  frozen  flow  is  the  lower  heat  flux  limit. 
We  note  that  the  discrepancies  with  the  results  obtained  assuming  a  reacting  boundary 
layer  are  rather  small  and  decrease  as  the  wall  temperature  increases.  Then,  for  wall 
temperatures  Tw  ~  1500/1  usually  encountered  in  catalycity  studies,  the  flow  is  nearly 
frozen  at  P  =  lOOOOPa.  Figure  34  also  shows  the  heat  flux  calculated  using  the  kinetic 
scheme  of  Dunn  and  Kang  and  Fick’s  law  for  diffusion  fluxes.  In  this  case,  the  influence 
of  the  modeling  of  the  diffusion  flux  appears  to  be  rather  small  and  the  use  of  Fick’s  law 
slightly  overestimates  the  calculated  wall  heat  flux. 


- N2 

- 02 

- NO 

- N 

-  0 


012345678 

11 

Figure  33:  Species  concentrations  in  the 
boundary  layer  for  Tw  =  1500/1  using 
the  phenomenological  model.  The  kinetic 
scheme  of  Dunn  and  Kang  is  used. 


Tw  [K] 

Figure  34:  Wall  heat  flux  calculated  for 
different  gas  kinetic  schemes  and  with  the 
kinetic  scheme  and  Dunn  and  Kang  and  the 
simple  Fick’s  law  for  diffusion  fluxes. 
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5.3.2  Low  pressure  test-case  :  P  =  2000 Pa 

We  consider  typical  boundary  edge  conditions  at  Pg  =  2000 Pa  :  Tg  =  5400 K  and  (3  = 
MOOOs^1  (Bourdon  et  al.  (2003)). 


Figure  35:  Wall  heat  flux  as  a  function 
of  the  wall  temperature  for  different  values 
of  the  average  7  coefficient.  The  kinetic 
scheme  of  Dunn  and  Kang  is  used. 


Figure  36:  same  as  Fig.  35  with  the  heat 
flux  calculated  using  the  phenomenological 
model  derived  by  Nasuti  et  al. 


Figure  37:  7  coefficients  of  the  phe¬ 

nomenological  model  derived  by  Nasuti  et 
al.  for  the  conditions  of  Fig.  36. 


Figure  38:  Different  values  of  0  coefficients 
for  the  same  condition  as  Fig.  37. 


Figure  35  presents  the  wall  heat  flux  as  a  function  of  the  wall  temperature  for  different 
values  of  7  and  Fig.  36  shows  the  flux  obtained  using  the  phenomenological  model  derived 
by  Nasuti  et  al.  As  for  the  high  pressure  case,  we  note  a  strong  nonlinear  dependence  of 
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the  heat  flux  with  the  recombination  probability  7.  The  heat  flux  for  a  fully  frozen  flow 
(i.e.  no  reactions  in  the  gas  and  at  the  wall)  is  also  shown.  The  difference  of  this  flux  with 
the  one  denoted  7  =  0,  puts  forward  (as  already  observed  for  the  high  pressure  case)  that 
for  P$  =  2000Pa  in  the  noncatalytic  case,  the  influence  of  the  chemistry  on  the  heat  flux 
decreases  as  the  wall  temperature  increases.  Figure  37  shows  that  7  coefficients  of  the 
different  species  calculated  with  the  phenomenological  model.  The  coefficients  7 00 1  Ino 
and  7 on  have  a  similar  evolution  with  temperature.  They  increase  up  to  Tw  ~  1500/1 
and  decrease  for  higher  temperatures.  For  300/1  <  Tw  <  2000/1 ,  ^/nn  increases  with  the 
temperature  and  for  Tw  >  1700/i ,  this  is  the  highest  recombination  coefficient.  We  note 
that  the  7 on  coefficient  is  much  higher  in  this  condtion  at  Ps  =  2000Pa  than  for  the 
previous  case  at  P5  =  lOOOOPa.  Figure  38  shows  the  surface  coverage  for  O ,  N  atoms 
and  the  OH  molecule  for  the  conditions  of  Fig.  37.  The  surface  coverage  by  the  OH 
molecule  decreases  as  the  wall  temperature  increases.  The  surface  coverage  by  O  atoms 
is  rather  high  for  300/1  <TW<  1500/1  and  then  decreases  rapidly  for  higher  temperature 
due  to  thermal  desorption.  The  N  atom  surface  coverage  remains  small  in  the  whole 
temperature  range  but  is  much  higher  for  Tw  <  1000/1  in  this  case  than  in  the  previous 
case  at  Ps  =  lOOOOPa  . 


Tw  [K] 

Figure  39:  Wall  heat  flux  calculated  for  different  gas  kinetic 
schemes  and  with  the  kinetic  scheme  and  Dunn  and  Kang  and  the 
simple  Fick’s  law  for  diffusion  fluxes. 


Figure  39  compares  the  wall  heat  flux  calculated  with  the  phenomenological  model  pro¬ 
posed  by  Nasuti  et  al.  using  different  kinetic  schemes  (Dunn  and  Kang  ,  Park  and  Gupta) 
and  for  a  frozen  flow.  The  heat  fluxes  calculated  using  the  kinetic  scheme  of  Dunn  and 
Kang  are  very  close  to  the  results  obtained  for  the  scheme  of  Gupta  for  the  whole  tem¬ 
perature  range.  The  heat  fluxes  calculated  using  Park’s  kinetic  scheme  for  Tw  >  1000/1 
are  slightly  higher  than  those  obtained  with  Dunn  and  Kang  scheme  but  the  discrepancy 
remains  small.  The  heat  flux  calculated  assuming  a  frozen  flow  is  the  lower  heat  flux 
limit.  We  note  that  the  discrepancy  with  the  results  obtained  assuming  a  reacting  bound¬ 
ary  layer  is  of  about  20%  for  Tw  =  300 K,  15%  for  Tw  =  1300/1  and  decreases  for  higher 
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temperatures.  We  note  that  the  influence  of  the  gas  chemistry  on  the  calculated  heat  flux 
appears  to  be  more  significant  for  the  boundary  layer  conditions  at  P$  =  2000Pa  than 
at  Ps  =  lOOOOPa.  Therefore,  for  wall  temperatures  Tw  ~  1500/1  usually  encountered  in 
catalycity  studies,  the  chemistry  of  the  gas  has  a  non  negligible  influence  at  Ps  =  2000Pa 
on  the  catalycity  coefficients  derived  for  the  material.  Figure  39  also  shows  the  heat  flux 
calculated  using  the  kinetic  scheme  of  Dunn  and  Kang  and  Fick’s  law  for  diffusion  fluxes. 
In  this  case,  the  influence  of  the  modeling  of  the  diffusion  flux  appears  to  be  rather  small 
and  the  use  of  Fick’s  law  slightly  overestimates  the  calculated  wall  heat  flux. 

Figure  40  shows  the  species  concentrations  for  Tw  =  1500/1  and  the  phenomenological 
model  proposed  by  Nasuti  et  al.  We  note  that  the  evolution  of  the  different  species  con¬ 
centrations  in  the  boundary  layer  is  very  different  than  the  one  obtained  for  lOOOOPa. 
Figure  41  shows  the  species  concentrations  for  the  same  boundary  conditions,  but  a  frozen 
chemistry.  We  note  that  a  15%  difference  on  the  wall  heat  flux  (Figure  39)  corresponds  to 
significant  discrepancies  on  the  different  species  evolutions  in  the  boundary  layer.  It  would 
be  very  interesting  to  compare  these  results  with  species  concentrations  measurements  in 
the  boundary  layer. 


Figure  40:  Species  concentrations  in  the 
boundary  layer  for  Tw  =  1500/1  using 
the  phenomenological  model.  The  kinetic 
scheme  of  Dunn  and  Kang  is  used. 


Figure  41:  Species  concentrations  in  the 
boundary  layer  for  Tw  =  1500/F  using  the 
phenomenological  model.  The  boundary 
layer  is  assumed  to  be  frozen. 


6  Conclusions 

Two  main  conclusions  can  be  derived  from  the  work  developed  in  the  previous  sections 
on  (1)  the  gas  phase  chemistry  and  on  (2)  the  interaction  between  a  gas  and  a  surface 
during  the  re-entry  of  a  spacecraft  into  the  upper  layers  of  the  earth’s  atmosphere. 

First,  the  elaboration  of  a  relevant  and  detailed  kinetic  scheme  for  a  given  gas  phase  is  a 
difficult  task  which  has  to  be  done  carefully,  since  numerous  species  and  elementary  pro¬ 
cesses  are  involved.  For  each  process  a  rate  coefficient  has  to  be  found  and  it  is  important 
to  try  to  find  the  most  accurate  theoretical  or  experimental  data.  The  chemistry  being 
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time-dependent  by  nature,  it  is  important  to  take  into  account  species  that  may  play  a 
role  in  the  time  evolution  even  if  their  concentration  at  the  stationary  state  is  negligible. 
Furthermore,  it  is  important  to  be  able  to  reduce  the  kinetic  scheme  after  having  carefully 
compared  the  results  of  the  simplified  scheme  with  the  detailed  scheme. 

We  have  used  this  procedure  to  derive  the  three-body  recombination  and  effective  ioniza¬ 
tion  rate  coefficients  of  atomic  oxygen  from  a  detailed  collisional  radiative  model.  The 
time-dependent  approach  used  in  this  work  has  clearly  put  forward  the  necessity  that 
the  system  be  in  a  quasi-steady  state  to  derive  meaningful  rate  coefficients.  In  addition, 
we  have  also  shown  that  the  value  of  these  rate  coefficients  depends  on  the  number  of 
excited  levels  considered.  Then,  to  elaborate  a  CR  model  for  air,  we  have  shown  that 
it  is  necessary  to  take  into  account  a  sufficient  number  of  excited  levels  in  order  to  cal¬ 
culate  correctly  the  ionization  or  recombination  as  well  as  all  processes  occurring  in  this 
type  of  medium.  The  results  obtained  put  forward  significant  discrepancies  with  common 
simplified  kinetic  schemes  used  for  re-entry  studies.  The  CR  scheme  developed  in  this 
work  can  be  used  to  derive  up-dated  and  new  simplified  kinetic  schemes  for  air  for  various 
applications.  For  the  validation  of  a  simplified  kinetic  scheme,  we  recommend  careful 
comparisons  of  the  time-dependent  results  obtained  with  the  simplified  scheme  and  with 
the  CR  model. 

Second,  the  stagnation  point  boundary  layer  study  has  shown  that  the  problem  of  the 
determination  of  the  catalycity  of  the  surface  interacting  with  a  re-entry  plasma  is  far 
from  being  solved.  The  surface  chemistry  plays  a  crucial  role  and  the  net  heat  flux  is 
highly  dependent  on  the  surface  temperature.  We  have  shown  that  the  chemistry  of  the 
gas  phase  plays  a  non  negligible  role  in  particular  for  P  <  lOOOOPa.  For  the  conditions 
studied  in  section  5  the  discrepancies  between  different  kinetic  schemes  remain  small,  as 
the  temperature  in  the  boundary  layer  remains  less  than  5400/F.  However,  the  time- 
dependent  CR  model  has  put  forward  that  these  discrepancies  may  be  significant  when 
the  temperature  in  the  boundary  is  high  enough  to  have  a  non  negligible  ionization  degree. 
In  this  work,  we  have  used  a  phenomenological  model  for  wall  catalycity.  We  have  shown 
that  this  type  of  model  gives  very  interesting  information  on  the  evolution  of  the  catalyc¬ 
ity  with  the  wall  temperature.  It  would  be  interesting  to  further  improve  the  modeling 
of  the  gas  surface  interaction  in  flow  and  boundary  layer  codes.  For  example  it  could  be 
very  interesting  to  take  into  account  the  dynamics  of  the  lattice  atoms  according  to  the 
semi-classical  collisional  method,  which  has  been  developed  and  successfully  applied  for 
the  description  of  oxygen  recombination  on  silica  surfaces  (Cacciatore  et  al.  (1999)). 
Finally,  to  validate  the  results  obtained,  a  comparison  with  experiments  is  necessary. 
However,  it  is  difficult  to  measure  the  concentrations  of  a  large  number  of  species  in  a 
flow  simple  enough  to  be  a  good  test-case  to  validate  accurately  a  nonequilibrium  chem¬ 
istry.  For  catalycity  studies,  up  to  now,  only  a  limited  number  of  parameters  have  been 
measured  in  the  stagnation  point  boundary  layer  and  at  the  wall.  It  would  be  interesting 
to  carry  out  further  measurements  of  the  evolutions  of  different  species  concentrations  in 
the  boundary  layer  to  better  validate  the  results. 
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